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NOMENCLATURE * 




m 

N 

n 

Pr 

P 

Pw 


2 

burning surface of propellant, cm 

2 

cross-sectional area of the port, cm 

2 

motor nozzle throat area, cm 

pre-exponential factor in burning rate law, ap n 
points at left boundary and right boundary, see 
Figs. 3-3 and 3-4. 
burning perimeter, cm 
discharge coefficient 

thrust coefficient accounting for losses 
thrust coefficient 
nozzle coefficient for thrust loss 
sonic velocity, cm/sec 

characteristic velocity in nozzle discharge relation- 
ship, cm/sec 

specific heat at constant pressure, cal/g-K 
hydraulic diameter of the port, 4Ap/b, cm 
total stored energy, cal/g 
thrust, newtons 

2 

friction coefficient, 2gt /pu 

2 w 

mass flow rate, pu, g/cm -sec 
conversion factor 

I . 

enthalpy of gases, cal/g 

local convective heat-transfer coefficient, cal/ 

2 

cm -sec-K 

inhomogeneous terms in the governing equations 
mechanical equivalent of heat 

3 

erosive burning constant in Eq. (2-20) / cm -K/cal 

Mach number 

mass flow rate, g/sec 

number of spacewise mesh points 

pressure exponent in the non-erosive burning rate law 
Prandtl number 
static pressure, atm 

wetted perimeter (after ignition, only the inert 
perimeter) , cm 


♦Appendices A and B provide further definitions of the 
symbols and units. 
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T 
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T ig 

T P i' T o 

T 

ps 

T 

ps,ig 


t 

u 

V 

W 

w 

X 

X 


seg 


slot 


e 


a 


N 


ex 


2 

rate of heat transfer, cal/cm -sec 
gas constant for the combustion gases 
universal gas constant 

burning rate of the solid propellant, including 
the erosive burning contribution, cm/sec 
burning rate under nonerosive conditions, cm/sec 
temperature (without subscript, static gas temper- 
ature) , K 

average film gas temperature, (T + T )/2, K 

ps 

adiabatic flame temperature of the solid propel- 
lant, K 

effective mean temperature of the igniter gas, K 

initial propellant temperature, K 

propellant surface temperature, K 

propellant surface temperature, at which propel- 

lant ignition occurs, K 
time, sec 

gas velocity, cm/sec 
volume, cm 3 

molecular weight of the combustion gases, g/g-mole 
width of opening to slot, cm 
axial distance from head end, cm 
effective distance from the leading edge of the 
corresponding segment, cm 

distance from head end of motor to the upstream 
edge of the segment, cm 

perpendicular distance from the propellant surface 
into the solid, cm 

2 

thermal diffusivity, cm /sec 

pre-exponential factor in erosive burning relation- 
ship 

nozzle half-angle at exit plane, deg 
erosive burning exponential factor 
ratio of specific heats 
nozzle exit expansion ratio 
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e = equivalent sand roughness, cm 

s 

8 = weighting parameter for implicit numerical scheme 

X = thermal conductivity, cal/cm-sec-K 

X N * nozzle divergence loss 

p = viscosity of the combustion gases, g/cm-sec (poise) 

3 

p = density (without subscript, gas density) , g/cm 

-1 

o = temperature sensitivity of burning rate, K 

r 

t = shear stress on the port wal) 

w 

<p - temperature gradient at propellant surface, K/cm 


Subscripts 
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am 

ch 
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gen 
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ig 

in 

int 
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mtr 
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out 

P 

pd 

pr,c 

ref 

sldv 

slot 

stag 


= values at points A 2 and A 2 , see Figs. 3-3 and 3-4. 
= actual 

* ambient 
= chamber 

= entrance to motor nozzle (motor chamber aft-end) 

= entrance section at head-end of motor 
= exit of nozzle 
= generated in slot 
= head-end 
= initial value 

* igniter 

= into slot 

* initial 
= measured 
= motor 

= index for mesh points in the axial direction 
= nozzle end 
= out of slot 
= pressure 

* latest predictor calculation 

= solid propellant (condensed phase) 

= reference conditions, Tq = 298K and p = 68.08 atm 
= per unit volume containing slot 
= slot 
= stagnation 
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I = 

II 

III 

1 , 2 , 3 , 4 = 


temperature 

velocity 

along right-running characteristic 
along left-running characteristic 
along particle-path line 
station numbers 


Superscripts 


j = index for mesh time increments 

— = mean value 
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1.0 INTRODUCTION 

This report describes the application of analytical 
techniques which are being developed to predict the ignition 
transients of segmented, solid propellant rocket boosters. 
Particular attention is given to the large, solid-rocket 
booster (SRB) which is being developed for the Space Shuttle 
(e.g., thrust > 10^ lb and length > 100 ft). The booster is 
referred to as a segmented motor (as opposed to a mono- 
lithic motor) since it consists of several large motor 
segments which are joined together. The development of 
large SRB's is accompanied by a number of questions pertain- 
ing to optimizing performance, improving reliability, and 
reducing the costs of qualification tests. During the 
formative phases of development, it is natural to employ 
analytical techniques that yield a more comprehensive 
understanding of complex interactions between the igniter 
gas flow, heat transfer to the propellant, flame spreading, 
developing flow field, and erosive burning. Indeed, because 
of the costs of manufacturing and testing each solid rocket 
booster, predictions and design recommendations based on 
comprehensive analytical models can play an important role 
in defining the SRB configurations. 

The essential elements that are incorporated into the 
analysis are: (a) techniques to account for the strong 

contribution of erosive burning, coupled to the chamber 
gas dynamics; (b) ability to consider significant pressure, 
temperature, and velocity spatial gradients, encountered in 
high loading density motors, and their variation with time 
during the starting transient, and (c) calculation of 
propellant surface heat-up to ignition and flame spreading, 
coupled to both temporal and spatial changes of flow param- 
eters in the chamber. Indeed, to account for the experi- 
mentally observed steep pressure rise and pressure over- 
shoots, the dominant influences of erosive burning and the 
spatial gradients of pressure, temperature, and velocity 



- 2 - 


must be accurately taken into account. Thus, using the tech- 
niques described in this report, analysts will be able to 
prescribe with greater confidence the upper limits of 
volumetric loading density and predict more accurately 
the pressure-thrust-time motor performance during the 
starting transient. 

While the emphasis of this report is on large SRB's, the 
methodology is not dependent on motor diameter and thus can 
be applied to a wide variety of high length-to-diameter 
ratio motors. 

The point of departure for this study was the analyti- 
cal solution and computer program described in Ref. 1. The 
analytical development in Ref. 1 is for monolithic motor 
configurations and is sufficiently complete to treat most 
high performance motors; however, the computer program 
described in Ref. 1 dealt with a rather narrow subset of the 
overall analytical formulation. Thus, the computer program 
described in Ref. 1 (while well suited for the special 
geometry, igniter, and burning conditions of the laboratory 
slab motor) is not suitable for analyzing high-performance, 
operational motors. Also, the Ref. 1 computer program was 
not prepared in modular form and, as a consequence, changes 
in geometry, burning rate laws, friction factor equations, 
convective heating equations, etc., are extremely difficult 
to perform. During this study these limitations were over- 
come by redoing the computer program and emphasizing the 
requirements that occur during the design of high perform- 
ance motors. 


TF£ 


gbi«inal 10 
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2.0 DESCRIPTION -OF MODEL 

2.1 Physical Situation Considered 

The emphasis of this study was on the type of 
segmented motor configuration and ignition events illustrated 
in Fig. 2-1. At the onset of igniter discharge, a flow 
pattern with pressure, temperature, and velocity distributions 
develops in the motor port. Subsonic flow with friction and 
heat transfer to the propellant at low pressure (usually 
between 1 and 3 atm for the low port-to-throat area ratios 
under investigation) is established in the chamber. During 
the early part of the transient, the motor nozzle flow is 
unchoked. The first phase of the starting transient is the 
induction interval (ignition lag) , which ends with the 
appearance of the first flame on the propellant surface 
(see Fig. 2-lb) . The ignition criterion (used in this study) 
is that a point on the propellant surface ignites when it 
attains some critical ignition temperature, typically between 
800 and 900 K. Thus, the process of flame spreading along 
the propellant surface, which starts at first ignition and 
ends when the entire surface is ignited, is described by the 
model of successive ignitions. Once started, the flame spread 
ing is accelerated by the increased heat flux due to the 
rapidly increasing mass flow originating from the already 
ignited propellant surface. As the flow encounters the 
separations between the segments, the boundary layer is 
disturbed producing local increases in the convective heating 
rates. Rapid chamber pressurization Usually begins with the 
onset of fast flame spreading. The last phase of the starting 
transient is the chamber filling period, which follows the 
completion of flame spreading. Prior to this time, the nozzle 
flow usually becomes choked. Erosive burning may take 
place during this period and thereafter, due to high gas 
velocities in the chamber. Also, the compression effects 
during the rapid pressurization will significantly increase 
the temperature of the chamber gases. Thus, the maximum 
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(a) Segmented rocket motor configuration 


SIGNIFICANT 



TIME 


(b) Significant ignition intervals. 

Fig. 2-1 Type of rocket motor and time intervals 
considered during study. 
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pressure may be much higher than the eventual equilibrium 
pressure. Significant pressure and temperature drops as 
well as velocity increases are established along the port. 
Following the pressure peak, a quasi-steady situation develops 
and pressure decreases toward an equilibrium value due to 
diminishing erosive burning, as the port area increases. 

2.2 Types of Ignition Transient Models 

Numerous studies have dealt with a particular funda- 
mental process taking place during the start of solid propel- 
lant motors or have focused on practical correlations to 
facilitate igniter design. A detailed review of previous 
studies is given in Ref. 1. The analytical models directed 
at the analysis of the overall starting transient, can be 
categorized into three major groups: (a) lumped chamber- 
parameter, or p(t) models (e.g., Refs. 2-6), (b) quasi-steady 

one-dimensional flow, or p(x) models (e.g.. Refs. 7 & 8) , and 
(c) temporal and spatial development of flow field (Ref. 1.). 

As shown on Fig. 2- 2a, the models of the first group 
assume uniform pressure and temperature distributions in the 
combustion chamber port. Accordingly, such models are 
incapable of considering those driving forces which control 
ignition and flame spreading rates. 

In the models of the second type (see Fig. 2-2b) , quasi- 
steady pressure, temperature, and velocity distributions 
along the port are considered at each instant of time during 
the transient. Flame spreading is often treated in ways 
unsuitable for predictive use, e.g., an experimentally 
determined function of time; partially or fully instantaneous; 
linear function of the burning rate; proceeding in a constant 
average rate; or calculated from an experimental pressure- 
time plot. 

Many high-performance solid propellant rocket motors have 
high volumetric loading densities, small port-to-throat area 
ratios, and, frequently, large length-to-diameter ratios. 

Such motors, referred to as HVT (High Velocity Transient) 
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ENERGY EQN. NOT USED 
T_,_ * T f + f [p (t) 3 


ch 

ENERGY EQN. 
T 


USED 


Ch 


= f (p(t) , V. . .) 

= f . . .) 


PROP 



INTERNAL GEOMETRY AND 
FLOW FIELD ARE NOT 
ACCOUNTED FOR ADEQUATELY. 



(a) p(t) model, chamber parameters are lumped. 



50% IGNITION 

I ^ 

10% IGNITION 


n 


CHAMBER VOLUME EFFECTS 
ARE NOT ACCOUNTED FOR 
ADEQUATELY . 


x 


SOLUTION IS A 
SERIES OF QUASI- 
STEADY STATE 
SOLUTIONS 


(b) p(x) model, quasi-steady 1-D flow. 


ENERGY, MOMENTUM 
& CONTINUITY EQ. 
APPLIED 


FLAME SPREADING 
RATE CALCULATED 




(c) p(x,t) model, transient flow field. 


Fig. 2-2 Several types of ignition models. 
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motors, are characterized by high internal gas velocities 
(Mach numbers > 0.3) , significant axial pressure and tempera- 
ture gradients, pressure over-shoots, and relatively short 
ignition transient times. The large solid boosters are 
properly classified as HVT motors. The starting transients 
of the HVT motors (and SRB's) are inadequately predicted and 
incompletely analyzed by either one of the above-mentioned 
classes of approaches. In particular, pressure, temperature, 
and velocity variations with time and axial position in the 
motor [p(t,x) model] must be considered (see Fig. 2-2c) . 
Interactions between processes, such as the developing flow 
field, igniter gas flow, convective heat transfer to the 
solid propellant, flame spreading, and erosive burning must 
be taken into account. 


2.3 Description of Model Used in This Stud^ 


The analytical model [referred to as the p(x,t) 
model] is capable of treating the following: 

General - 

‘Goes from onset of igniter flow through the ignition 
transient to steady state operation. 

Igniter - 

‘Igniter discharge is an arbitrary function of time. 

‘Pyrogen igniter located at head-end. 

Main Chamber - either monolithic or segmented configurations: 

‘Axial gradients of pressure, temperature, and velocity. 

‘Nozzle flow either choked or unchoked. 

‘Axial variation of port area. 

‘Wall friction and heating losses. 

‘Propellant surface temperature and gradient varies 
in axial direction. 

‘Flame spreading is governed by coupling between main 
chamber flow field, convective heating rates, and 


propellant temperature profile. 
‘Erosive burning is accounted for. 
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The key to quantitative P max predictions for SRB 
designs is expected to be an accurate determination of the 
erosive burning during chamber pressurization. 

Even though the physical model and mathematical formula- 
tion given in Ref. 1 were rather comprehensive with respect to 
high performance rocket motors of monolithic configuration, 
the Ref. 1 computer program was very specialized and suitable 
only for analyzing the laboratory slab motor. Thus, the 
program was not sufficiently general to handle conventional 
rocket motor designs. Accordingly, during this study the 
computer program was reorganized, improved greatly, and 
expanded to give it the capability to analyze high performance 
rocket motors. Thus, the reader may refer to Ref. 1 for the 
details of the equation development for monolithic motors. 

Of course, the extensions required to consider the slots 
between the segments are documented in this report. The 
duplication between Ref. 1 and this report has been kept to 
a minimum. Furthermore, the Ref. 1 nomenclature has been 
retained. 

2.4 Basic Assumptions 

The following basic assumptions are used in the 
p(x,t) model: 

(a) All chemical reactions occur on the propellant 
surface in a combustion zone which is so thin 
that it can be considered as a plane. The combus- 
tion products enter the main stream in the port 
with zero axial momentum. 

(b) The flow in the main chamber port is one-dimensional. 
This assumption is a good approximation to the real 
situation in view of the highly turbulent flows 
encountered in HVT motors. Changes of flow proper- 
ties across the boundary layer are considered in 

the expressions for heat transfer and friction co- 
efficients used in the analysis. It is further 
assumed that the slots are designed so that their 
interactions do not affect the one-dimensionality 
of the flow. 
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(c) Rate processes at the propellant surface are quasi- 
steady in the sense that their characteristic times 
are short compared to that of the pressure transient. 

(d) The propellant combustion products and the pyrogen 
igniter gas have the same values of c^, W, and y. 

(e) The gases flowing in the port obey the perfect 
gas law. 

2 . 5 Governing Equations 

The mathematical formulations of the approach to 
the above described problem consist of the following: 

(a) Mass, momentum, and energy conservation equations 
in unsteady, quasi-one-dimensional form for the 
gas phase; 

(b) Equation of state for the gas flowing in the motor; 

(c) Proper initial conditions at the start of the 
transient (onset of igniter flow) ; 

(d) Two boundary conditions at the fore-end of the 
propellant section, obtained from a pair of 
ordinary differential equations, which describe 
the rate of change of pressure and temperature 
in the entrance section; 

(e) A third boundary condition, which describes the gas 
velocity at the entrance to the motor nozzle, for 
either choked or unchoked flow; 

(f) Semi-empirical correlations for the convective 
heat-transfer and friction coefficients for the 
highly-turbulent flow in the port; 

(g) Burning rate law for the solid propellant, including 
the effects of initial temperature, pressure, and 
velocity over the surface (erosive burning) . (The 
burning rate equations are in subroutine form so 
that they can be easily modified.) 

(h) A solid phase heat-up equation for determination 
of the propellant surface temperature during the 
induction interval coupled to an ignition criterion 
for the propellant; and 
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(i) Thrust calculations based on nozzle-end stagnation 


pressure. 

First, the conservation equations along a section of the 
port that does not contain slots will be summarized. For a 
more detailed development, the reader is referred to Ref. 1. 
The mass conservation equation for the control volume in Fig. 


2-3, written for a unit volume, is 



+ 


9 (pu) . 
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The momentum conservation equation is 
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The energy conservation equation, written in terms of the 
, total stored (internal and kinetic) energy per unit mass, E, 
is 


9 

Jt 


1 


3 

• 

[ pA p E J 

+ 


P V E , 


RATE OF ACCUMULATION OF RATE OF SPATIAL INCREASE 
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NET RATE OF 
ENERGY INPUT 
BY HEAT 
CONDUCTION 


RATE OF WORK 
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RATE OF WORK DONE 
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VISCOUS FORCES 


V rh fPw “ q sP W 


RATE OF ENTHALPY 
ADDED BY BURNING 
OF SOLID PROPELLANT 


RATE OF HEAT 
LOST TO THE 
SURROUNDINGS 


(2-3) 


The above conservation equations constitute a set of 
inhomogeneous, nonlinear, partial differential equations (PDE) 
After nondimensionalization and an order of magnitude analysis 
the effects of the following terms were neglected (indicated 
in Eqs . 2-2 and 2-3): (a) forces between gas molecules 

due to the viscous stress in the axial direction, created 
by the axial velocity gradient; (b) viscous dissipation 
and rate of work done by the internal viscous forces, and 
(c) axial heat conduction between gas molecules. The wall 
friction at a station in the propellant section is neglected 
after ignition occurs and appreciable blowing starts. 

The present formulation is more general than the Ref. 

1 formulation since it accounts for the rate of mass accumu- 
lation in the free volume created by the propellant surface 
regression. 


2.6 Provisions for Circumferential Slots in Segmented 
Motors 

As shown in Fig. 2-1, segmented motors consist of 
a series of sections which when joined together form circum- 
ferential slots which contain an appreciable volume and 
burning surface area. The interactions produced by the slots 
are accounted for by solving a set of ordinary differential 
equations for each slot and using those solutions as boundary 
conditions for the main stream flow equations. 
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Two flow interaction conditions between the main chamber 
and the slots are considered: (1) when the pressure gradi- 

ent between the main chamber and slot is significant (such 
as during the passage of the first strong pressure wave) , 
and (2) when the pressures in chamber and slot are nearly 
equal. In principle, the first condition could be used 
throughout the ignition transient but to do so greatly 
increases the computer time without increasing the accuracy 
of the solution. 

2.6.1 Pressure Differential Between Main Chamber 
and Slot is Significant 

Whenever the pressure differential between the 
main chamber and the slots is significant, the coupling of 
the fluid flow processes between the several circumferential 
slots and the main chamber is accounted for by simultaneously 
solving the continuity, energy and momentum equations for 
flow into (or out of) each slot and coupling the result with 
the complete PDE solution for the main chamber flow. 

In Fig. 2-4, the control volume of the slot is shown 
interacting with the main stream flow. The flow in (or out 
of) the slot is described by the following equations. 

When the main chamber pressure is greater than the slot 
pressure, the flow into the slot is 


m. = C^bw . ,p ' 

in D slot H ch 


2 9 Jc p T ch 


(^slot^chj 


]2/y 


p slot /p ch 


(Y+1)/y 


1/2 


(2-4) 


and 


m , = 0 
out 


Conversely, as the slot begins to generate mass and the pres- 
sure in the slot exceeds the pressure in the main chamber, 
flow is out of the slot. 
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MAIN STREAM FLOW 



T ch' P ch ,P ch 


m 

i , out 

^ / 

m 


BURNING AREA, 

6A, 1 ^ 

b, slot 


slot' p slot' p slot if 


AREA CONNECTING 

SLOT AND PORT 

w , , b 
slot 



INDICATES 

INHIBITED 

PROPELLANT 


slot" 


\ \ \ 


Fig. 2-4 Longitudinal section through slot showing 
nomenclature used in analysis of slot. 
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m out ^D^slot^slot 


2gJc T , . 
3 p slot 


^ch^slot 


2/y 


p ch^slot 


(Y+l)/Y 


1/2 


(2-5) 


and m. = 0 
m 


After ignition, the mass generated in the slot is 

m gen ~ r b p pr^b, slot 


( 2 - 6 ) 


The rate of increase of slot volume is 

dV slot/ dt = r<A b,slot - w slot b) 

The rate of mass increase in the slot is 

d ^ V slot p slot^ dt - ” m out + m in + m gen 


(2-7) 


(2-8) 


where either m. or in . is zero, 
in out 


In terms of the sign convention for the main chamber equation 


m i . = in . - m. 
slot out in 


(2-9) 


The pressurization rate of the gases in the slot is 
dp slot /dt = 


Jc [m . T - . - m. T , - m T-. + p . . (d(V). , ../ dt)/R] 
p out slot in ch gen f ^slot slot ' 1 


V slot (1 - C P J/R) 


( 2 - 10 ) 


Gas temperature in the slot is calculated from the perfect 
gas law. 

2.6.2 Pressures in Slot and Chamber Nearly Equal 

After the pressures in the main chamber and slot 
become nearly equal, the coupling of the fluid flow processes 
between the slot and main chamber is accounted for by simul- 
taneously solving the energy and continuity equations for the 
flow into or out of each slot. 
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Combining the energy and continuity equations yields 
an expression for the rate of temperature change of the slot, 


dT 


slot 

dt 


rt'I 

■ 

taJslot 

m 

gen 


yT,. - T , .1 + m. 

T f slotj m 


YT ch T slot 


dV 


" Ia out (Y ~ 1)T slot " (Y “ 1) p slot T slot'dt 


slot 


( 2 - 11 ) 

Differentiating the equation of state produces an expression 
for rate of mass increase in the slot 


d ‘ pV »slot 

dt 


E v_ 

RT 2 


dT 


slot 


dt 


slot 



slot 



lH slot dt 


( 2 - 12 ) 


In terms of the sign convention for the main chamber 
equations 


m. 


d ' PV) 3lOt 


+ m 


slot dt "gen 

When m s i ot > 0 • net flow into main stream. 


( 2 - 13 ) 


When m s i 0 t £ net flow out of main stream. 

The above set of equations which specify p, pV, V, 
and T for a slot were solved using a fourth order Runge- 
Kutta Method at each At. 


2.7 Coupling Slot Flow to Main Chamber Flow 

The elements considered in the main chamber continu- 
ity, momentum, and energy equations are illustrated in Figs. 
2-5a, 2-5b, and 2-5c respectively. An order of magnitude 
analysis revealed that for pressures and accelerations being 
considered, viscous normal stresses and axial heat conduction 
in the gas produced relatively small effects and could be 
neglected with no loss in accuracy. The derivation of the 
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Fig. 2- 5a Elements considered by 

continuity equation. 
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Fig. 2-5b Elements considered in 
momentum equation. 
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p r c T^bAx ’ 
pr d p t 



Fig. 2-5c Elements considered in energy equation. 
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interaction of the slot flow with the main chamber flow yields 
the following governing equations for those increments which 
contain slots. 


3u 3u , RT 9£ _ ~ p pr r b b RTu _ m sldv RTu 
at U 3x Sx Ap p p 


fp u 2 
w 

2A 

P 


(2-14) 


3T + U 3T + ( 1)t 3u 

at 3x ; 3x 


, . , „ 3A r, b „ 

A 9x r M pr A d 

p P ‘ 


m T u 

T - — + — 

f y 2gJc 


Pi 


JTra s lt av^ 

+ (Y - 1) „ sx ih 


u 


h , Slot 

[ slot 2gJ 


RT 


(Y-l)JJ 


- (Y 


1 ) 


p T 




, 3 ) 
fpu 

2gRT, 


(2-15) 


l£ 

at 


+ u|£ + 

3x 


YP 


, 9u 

ax 


A 


a a 

— 1 
ax 


p r, b 

• 

f 'l 

2’ 

+ / v _ i ) P r b 

i 

T - T P 

+ iL_ 


C u 

p 

lf Yf' 

2? 

p 


P r J 



where m s ^ v is mass flux from slot per unit volume of 
increment containing slot, e.g., m ,i ot /^p &x . 


fpu 3 
£ ~ 2gRT j 

(2-16) 


Equations (2-1 4) , (2-15) , and (2-16) for the segmented 
configuration are the counterparts of Eqs. (4), (5), and (6) 
respectively in Ref. 1 for monolithic configurations. The 
primary difference in the two sets of equations is the 
addition of inhomogeneous terms to each equation when slots 
are considered. 
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2,8 Empirical Correlations for the Heat-Transfer and 
Friction - Coefficients 


The expression for the local convective heat-transfer 

coefficient h is deduced from the conventional Dittus- 
c 

Boelter correlation for turbulent flow in pipes. Entrance 

effects, observed by many investigators (e.g., Refs. 9-13) 

are accounted for by including a power function of the length 
1 12 

and diameter. ' Variation of the gas physical properties 
across the boundary layer is considered by evaluating the 
properties at an average film temperature T a f The expres- 
sion for h c used in this study, is 



1.56 x 10" 3 Pr~°* 6 c 

P 


pul 0.8 0.1 -0.67 

Ir 1 W r af 


x seg d hJ 


- 0.1 


(2-17) 


where x is the effective distance from the leading edge 
seg 

of the segment being considered. To avoid complications that 

arise when unrealistically small values of characteristic 

length are used in convective heating relationships, the 

minimum value of x se g is considered to be one-half of a 

distance increment, Ax, as used in the numerical solution. 

After ignition, convective heating to the propellant surface 

is terminated. However, the calculation of h is continued 

c 

for use in the erosive burning calculations. 

The friction coefficient, f, is deduced from Colebrook's 

expression for turbulent flow in pipes with roughness. 

14 

Observed entrance effects and variation of fluid properties 
across the boundary layer are taken into account. The friction 
factor is expressed in the following form: 


f = 


0.449 (^ h /x seq) 


0.1 


/ j •) ___ - 0 . 5 m 1 • 65 x 0 . 05 

V ln | £ s /d h + T af (d h /x seq ) 

3,7 10 6 pud h f 0,5 


(2-18) 
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2.9 Burning Rate Law 

Two versions of the Lenoir-Robillard burning rate 

15 X 6 

law ' can be used to account for erosive burning. The 
first is the conventional form 


r - r 0 ♦ V=°- V ' 2 ~P 


- I B e rp pr /G 


(2-19) 


h is the hydraulic diameter, 4A p /b, at the point 


where d 

being considered. Equation (2-19) accounts for the convec- 
tive effects by incorporating the important functional 

dependencies on free stream flow and geometry, i.e., 

„ 0 . 8 X . 0.2 

G /d h . 

The second is a form which is particularly well suited 
for analyzing segmented motors. 


r = r„+ k e h c (d lf x,L slot ,u,...)exp 


- 8 rp /G 
1 e K pr' 


( 2 - 20 ) 


where h c is any convective heat transfer coefficient which 

accounts for the intricacies of the solid rocket geometry 

and flow field, in particular, those produced by the segments. 

In principle, k and a can be estimated from 

e is e 

theoretical considerations, but in practice they are 
usually deduced from experimental data. Values of 3 e are 
always determined from experimental correlations and general- 
ly are in the range between 75 and 200. In all cases, it is 
necessary that the k and 3. pairs and the a and 3 

c c c c 

pairs be consistent. 


2.10 Determination of the Propellant Surface Temperature 

The heat equation for an unignited propellant grain 
at each of the series of axial locations (e.g., 20 to 29) is 


3T (t,y) 3 2 T (t,y) 

"at — v — 


ay 


with the initial condition 


( 2 - 21 ) 


Tpr ( 0.y> = T pi 


( 2 - 22 ) 
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and with the following boundary conditions: 

T pr (t ’”> = T pi < 2 - 23a > 

3T pr (t,0)/3y = -(h c / X pr ) (T ' T pr ) (2-23b> 

Equation (2-21) along with its boundary conditions is solved 

17 

using Goodman's integral method. 

For purposes of this study, ignition of a particular 
section of the propellant surface is accomplished when a 
prescribed critical ignition temperature is attained. This 
simple ignition criterion is satisfactory whenever high 
convective heating rates assure rapid ignition. However, 
if ignition failures (i.e., hang fires and misfires) are 
to be studied, a more comprehensive ignition model may be 
required . 


2.11 Head-end and Nozzle-end Boundary Conditions 


The three physical boundary conditions are as 
follows: two boundary conditions at the entrance to the 

propellant section (denoted by x ), and one boundary con- 

■LT 

dition at the entrance to the motor nozzle (denoted by x„) . 

ill 

The gas properties in the entrance section are assumed 
to be uniform. Heat-transfer and shock-pattern losses taking 
place at the exit of the igniter may lower the effective 
mean temperature of the igniter gas» T^ . From the contin- 
uity and energy conservation equations for the entrance 
section, the rate of change of pressure and temperature at 
the entrance to the propellant section are derived: 


dp 

r es 

dt 


es 


yRT. m. (t) - 
1 lg ig 


YA„ oc P ae (t)u (t) 
p,es es es 


- 1|A p,es p es tt)u es (t) 

2gRT es (t> 


(2-24) 
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dT 

es 

dt 


1 Ks<t)n. ia (t) 
n — P it) — 

es *es 


|> T ig - T es (t) ] - 


^ - 1)A p,es u es (t) 


T es (t) + 


u e s (t > 

2gR 


(2-25) 


The igniter mass flow rate m^ is a prescribed function 
of time. 

Equations (2-24) and (2-25) form a pair of coupled 

ordinary differential equations for the unknowns p and 

es 

T . Their solution provides two boundary conditions for 
es 

the governing Eqs. (2-14 to 2-16). The two initial condi- 
tions for this pair of equations are 


p es (0) = p es,i and T es {0) 


T 

es,i 


(2-26) 


The boundary condition at the aft end of the motor 
chamber is derived assuming isentropic flow between the 
entrance to the motor nozzle and the nozzle exit. Using 
continuity and energy conservation relations between these 
two sections, the following equation applies: 


U E (t)2 - f^ RT E (t) 


“ ^P ex (t) /p E (t) ] ^ ^ (2-27) 


(A /A eJ{ ) Ip ex (t) /P E ( fc ) ] 


When the motor nozzle is unchoked (the usual case during 
the induction intervals, when very low port-to-throat area 
ratios are used) the pressure at the nozzle exit, p in 

GX 

Eq. (2-27), equals the ambient pressure. For the case of 
a choked motor nozzle Eq. (2-27) becomes the following 
implicit relation: 


u E (t) 2 


ygRT E (t) 

vv i L ' 


2 

TfT 


l + 



(Y+1)/(Y-1) 
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3 ♦ 0 NUMERICAL SOLUTION 

This section extends the discussion of the numerical 
solution in Ref. 1 by describing the organization of the 
revised computer program and by describing the methods used 
to account for the slots. The overall numerical solution is 
summarized in this section. However, the reader is referred 
to Ref. 1 for a more general discussion of the specific 
numerical techniques which are used. 

3.1 Provisions for Circumferential Slots in Segmented 
Motors ~ ’ ~ 

As shown in Fig. 3-1, segmented motors consist of 
a series of sections which when joined together form circum- 

i- 

ferential slots which contain an appreciable volume and 
burning surface area. The interactions between the main 
chamber flow and the slots are accounted for by coupling the 
equations for the main chamber with the ordinary differential 
equations for each slot. During the development of the 
numerical solution, certain conventions regarding geometry 
were adopted: 

( 1 ) 


( 2 ) 

( 3 ) 


Segmented motors are divided into equal axial 
distance increments (Ax) regardless of whether an 
individual increment includes a circumferential 
slot. Circumferential slots provide an additive 
volume (6V g ^ ot ) and an additive burning surface 
area (SA^ slot ) but do not alter the port area, 
i.e., the mainstream flow does not pass through 
the slot. 


Port area A p is the part of the flow passage 


above the slot and A = A *. + brAt where 

p,t+At p,t 

b is the burning perimeter along the port. 

Burning surface area of an increment which includes 
a slot is 
the value 


AA, . , = bAx + 6A, , . . Note that 

b,slot b,slot 

6Afo slot must be input as the burning 


area in the slot less bw 


slot* 


The burning surface 



PORT AREA BURNING SURFACE AREA 


HEADEND 

SEGMENT 

WITH 

COMPLEX 

GEOMETRY 


COMMON 
SEGMENT 
NO. 1 


COMMON NOZZLE END 
SEGMENT SEGMENT 
NO. 2 


x 


E 





Fig. 3-1 Geometrical description of segmented rocket motor 

showing lumped parameter method of treating volumes 
and surface areas associated with slots. 


REPRODUCIBILITY OF TR 
jsRISfNAL PAGE IS POOR 
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in the slot may be quite small since the circumfer- 
ential area in the slot may be partially inhibited 
to prevent burning. AA^ slot is considered to be 
constant during the ignition interval. 

(4) The total volume of an increment which includes a 

slot is AV = A Ax + 6V . . . However, only A Ax 
p slot J p 

is associated with the main stream. During each 
At, the volume of the slot increases by rAt(6A^ s ^ ot 


- bw slot> • 

(5) The pre-ignition heat loss and friction are affected 
by the slots since the slots disturb the boundary 
layer and accumulate gases during pressurization. 

(6) Ignition of the exposed propellant in the slots is 
controlled by the flame spreading along the main 
port. Since in the present SRB design the burning 
surface area in the slots is a small fraction of 
the total burning surface area, the burning surface 
area in a slot is assumed to ignite at the same 
time as the axial increment that contains it. 

(7) The passage area into a slot is bw s p 0 t anci i n “ 
creases by 2rAt during each time interval. The 
mass from the slots is injected norma 1“ to -the - 
centerline. 

The tabular inputs required to describe segmented rocket 
motor geometry are: 

For both monolithic and segmented motors: 


x 


P 


' w 


Distance from head end of motor, i.e., 
x = 0 at head end of motor. 

Port area of main channel. 

Perimeter of propellant surface in main 
channel . 

Total perimeter of propellant and inert 
surface (i.e., wetted perimeter before 
ignition) . 
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For segmented motors, when slot begins in interval from 
x to x + Ax: 


‘slot 


*b, slot 


Distance to the upstream edge of the slot 
from the head end of the motor, i.e., x = 0, 
Burning surface area in slot. (The lack 
of surface area along the main port because 
the slot width is accounted for in this 


input value, i.e., subtract bw 


from 


slot 


w_ 


slot 

the actual burning surface area in the slot.) 
Volume (region of gas not in the flow field) 
of slot. 

Width of slot 


are the same for both 


slot 

The geometry inputs x r , x_ and 

P h La 

monolithic and segmented motors. Figure 3-2 illustrates the 
conventions used in defining the geometry. Note that the 

head-end region between x = 0 and x = x is a flow passage 

P 


that does not contain propellant and that the last 
a flow passage that does not contain propellant. 


Ax/2 is 


3.2 Implementation Scheme 

Numerical mathematical techniques were implemented 
for the simultaneous integration of the three governing 
equations, [Eqs. (2-14) , (2-15) , and (2-16)], the equations for 
the entrance section [Eqs. (2-24) and (2-25)] and exit section 
[Eqs. (2-27) or (2-28)], and the equation describing rate of 
change of the propellant surface temperature, [Eq. (III-50) 
in Ref. 1]. The numerical steps were organized in modular 
form into a program for a large capacity digital computer. 

The set of governing equations is totally hyperbolic in 
nature. ' All three eigenvalues of its characteristic 
equation are distinct and real. A generalized implicit 

20 21 

scheme based on central differences in spacewise derivatives ' 
was chosen to solve numerically the governing equations. Let 
the net of points in the t-x plane be given by t = jAt 
and x = nAx, where j = 0,1,2..., and n = 0,1,2... . The 
mesh size of the net is determined by At and Ax. Then the 
derivatives of pressure, for instance, are expressed in the 
following difference form: 
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FOR MONOLITHIC MOTORS - 

FROM INPUT TABLES OF PORT AREA AND PROPELLANT PERIMETER VS 
DISTANCE PROGRAM INTERPOLATES FOR b AND A AT EACH STATION 
1 THROUGH NDELX . P 

FOR SEGMENTED MOTORS - 
PORT AREA, PERIMETERS ARE INPUT FOR EACH Ax. 

x E ~ x G = Ax/2 



n=0 1 2 3 4 5 NMAX 


BPO BP (1 ) BP (2) BP (3) BP (4) BP (5) BP (NMAX) 

APO AP (1) AP (2) AP (3) AP (4) AP(5) AP (NMAX) 

PROPELLANT AREA: 

NDELX- 1 x 

A ^ 

A, = ~BPO + l AxBP(I) % bdx 

D 1 1 x 

1 P 

CHAMBER VOLUME: 

NDELX- 1 f X E 

V = x APO + f^APO + l AxAP(I) + I^AP (NDELX) % A dx 
P * ■, z p 

X J 0 

Fig. 3-2 Accounting for surface area and volume. 
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ft. 


3+1 j 

P - P 
n f_n 

At 


i£ 

3x 


fl Lj+! _ n j+l' 

(^n+1 ^n-1. 


n 


+ (1 


2Ax 


0) |p.^ 


f (3-1) 


n+1 


Pn-1 


) 


where * p(jAt, nAx) . 

The weighting parameter 6 is a real constant, lying 

in the interval 0=0=1. For the implicit formulation no 

20 

stability restriction exists if 0.5 = 0 = 1. Usually a 
value of 0 = 0.6 was used in the calculations. 


3.3 Quasi-linearization and Predictor-Corrector Calcula - 
tions of the Non-linear Terms in the Governing 
Equations 


The governing equations may be written in the 
following matrix form 


u 


(F 

u, u 

F. m 
U, T 

F 

u,p 


u 




T 

+ 

F 

T, u 

F 

T,T 

F 

T,p 

3 

3x 

T 



(3-2) 

P 

J 


F 

P/U 

F m 
p,T 

F 

P e P J 


P 

. 


I 

PJ 



where F[u(t,x) ,T(t,x) ,p(t,x) ] are the functional coefficients 
of the partial spacewise derivatives, and I [u (t,x) ,T (t,x) , 
p(t,x)] are the corresponding inhomogeneous terms in Eqs. 

(2-14) , (2-15) , and (2-16) . These latter terms account for 
the interactions with the slots. 

To obtain a system of linear algebraic difference equa- 
tions which can be solved simultaneously by matrix methods, 
the nonlinear coefficients F(u,T,p) and inhomogeneous terms 
I(u,T,p) are linearized in a way described in the following 
paragraphs . 
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The inhomogeneous terms were quasi-linearized as follows: 

n 


22 


p+e % I 


+ 0 

\ j+1 j] 

u J - U J I 

! 31 

i — 

n 

[ n n nj 


' n nj 

1 3u 


+ 0 


0+1 


n 


- T- 


1 31 

3 + 2 + J 

j 8T 

+ 6 
n 


j+l 


n 


- P 


n 


31 


3p 


j+ ! 


n 


(3-3) 


where the partial derivatives, which are actually algebraic 
terms, are evaluated with the flow properties at (j + 6/2, n) , 


such as, for instance, 


J+6/2 


n 


= p J + 0 


J+1 


n 


- p 3 


/ 2 . 


A predictor-corrector technique was applied to the non- 
linear coefficients F(u,T,p) and the partial derivatives in 
Eq. (3-3) . According to this method, a predictor calculation 
is first made, in which previous-time flow properties are 
used to evaluate all the coefficients F and partial 
derivatives in the quasi-linearized inhomogeneous terms. 

Then, a corrector calculation is made in which the coeffi- 
cients F and the partial derivatives are evaluated with 
weight-averaged flow properties. The iterations may continue 
until sufficient covergence is achieved. 

After substitution of the aforementioned finite-differ- 
ence representations into the governing equations, the 
velocity variation equation, Eq. (2-14) , becomes 


u j+1 - u j 

n n + U 3 

0 

j+1 j+1 

u n+l u n-l 

+ (1 - 6) 

-j j 

u 3 - u J . 

[ n+1 n-lj 

- 

At u n 

L 2Sx J 


T 

+ gR- 
P 


’0 

j+1 j+1 

[ P n+1 • P n~lj 

+ (1 - 0) 

f • . 

p n + l - p n-l. 


2Ax 


n 


u 3 ,T 3 ,p 3 , 
n n 1 nj 


+ 9 


( • _i_i . 1 3 1 

l u ] + l . u 3 “ 

n nj 3u 


3+ 2 
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+ 0 


t 3+1 _ T 3 
n n 
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3T 


j+- 
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+ 9 (p 3+ ^ - p 3 
' r n *nj 
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dp 


3+ ! 
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(3-4) 
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The temperature and pressure variation equations, Eqs. 
(2-15) and (2-16) are treated in an analogous manner. 

The barred parameters are weight-averaged quantities 
as follows, e.g., 


u^ = 0 
n 


u 


j+1 


n 


. + (1 - 0)u^ 
pd n 


(3-5) 


whereas the flow parameters in the partial derivative form 
in Eq. (3-4) are 


0 


j-Hy 9 
u n 2 = 2 


u 


j+1 


n 


pd 


1 - i 
x 2 


u- 


n 


(3-6) 


The parameters ( ^ ) ^ are the quantities calculated 

from the latest predictor iteration. In the first step 
predictor calculation they have the values of the previous 
time step. 

Negligibly small differences were found between final 
solutions obtained by single-step and multi-step predictor 
calculations. Therefore, to reduce computing time, single- 
step predictor-corrector calculations are used in the actual 
computations. The final solution at every time step is 
obtained implicitly by a matrix method discussed in 
Section 3.6. 


3.4 Extraneous Boundary Conditions 

The use of central-difference formulation, as 

described in Eq. (3-1) , for all spacewise derivatives in the 

governing equations, requires six boundary conditions for 

the solution. In other words, the set of finite-difference 

19 

equations corresponds to a hyperbolic system of 6th order. 
Therefore, three boundary conditions, in addition to the 
physical ones (defined by Eqs. (2-24), (2-25), and (2-27) 
or (2-28)] are needed. These so-called extraneous bound- 
ary conditions for the hyperbolic system of governing 
equations studied are derived from the compatibility 
relations at the boundaries. The latter are obtained 

by solving the governing equations by the method of char- 

. .. 23 

acteristics. 
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The compatibility relations along the right-running 
and left-running Mach lines (defined by dx/dt = u ± c, 
respectively) , in terms of p-u characteristics, are 





(3-7) 


The relation along the particle-path line (defined by 
dx/dt = u) , in terms of T-p characteristics, is 


dTl = Y-l T fdp 

, dt Jiii y p l dt J II;[ 



Y-l h 
Y P P 


(3-8) 


In relations (3-7) and (3-8), I u , I T , and I p are 
the inhomogeneous terms in the governing equations [see 
Eq. (3-2)]. 

For a subsonic flow toward the nozzle in the motor which 
is the case during most of the transient in the physical 
model,, described in Section 2, the left-running characteristic 
line is considered at the left boundary, whereas the right- 
running characteristic line together with the particle-path 
line are considered at the right boundary of the region of 
numerical computation. Taking into account the conditions 
at the right boundary, which is the aft-end of the motor 
(uniform perimeter, no area change, and no mass addition) , the 
compatibility relation along the right-running characteristic 
is the same as given in Ref. 1 as Eq. (IV-12) . The compati- 

bility relation along the left-running characteristic is the 
same as Eq. (IV-13) in Ref. 1. The compatibility relation 
along the particle-path line, at the right boundary, is the 
same as Eq. (IV-14) in Ref. 1. 


3.5 Treatment of the boundary Conditions 

The compatibility relation along the left-running 
characteristic line, Eqs. (2-24) and (2-25), form a 
closed system of ordinary differential equations to 
determine the gas velocity, temperature, and pressure at the 
entrance to the propellant section. The system, formed by 
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the other two compatibility relations, Eqs. (2-27) or 
(2-28) , (unchoked or choked flow, respectively) deter- 
mines the flow parameters at the chamber aft-end. The 
equations in each system are expressed in a finite-difference 
form and quasi-linearized in an implicit form, using a pre- 
dictor-corrector iteration calculation in the same manner 
which is utilized for the governing equations. As a result, 
a set of three linear inhomogeneous algebraic equations is 
obtained for every- boundary . Each set of equations is solved 
by 3 x 3 matrix, using Kramer's rule, to yield the gas 
velocity, temperature, and pressure at the boundary for the 
next time step calculations. The characteristic line 
segments near the left and right boundary of the calculated 
domain are shown in Figs. 3-3 and 3-4, respectively. These 
figures also show the numerical calculation grid in an 
x, t-diagram. 

In Fig. 3-3, the line sent out from A~ to the left- 

i+1 ^ 

boundary point B^(xQ,t J ) with a slope (dx/dt).^ = u - c 
is the left-running characteristic. The process of deter- 
mination of the boundary values is carried out in the 
following steps: 

(1) With the calculation completed for the time step 
jAt, all properties are determined at the point 

by linear interpolation between the boundary 
(x^t- 1 ) . The position of A;j is taken from the 
previous-time-step calculations, as described 
below. 

(2) The left-boundary system of equations [Eqs. (2- 

24), (2-25), and (IV-13) in Ref. 1] is solved 

i+1 

for u, T, and p at B^(xQ,t J ). The flow 
parameters, present in the coefficients and inhomo- 
geneous terms of these equations, are iterated 
with the solutions until sufficient convergence 
is reached. 
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(3) The corrected location of A^ at the current time 
is found from the slope of the left-running char- 
acteristic (solid line - B„ in Fig. 3-3), i.e.. 


j _ At 

X A 2 X 0 “ 2 


j+1 , j 

U x U A 
0 2 J 


c j+1 + <4 

x 0 fl 2 


I£ I x a 2 " X ? 2 I /X ; 

are repeated with (x 
subscript 


(3-9) 

> 0.002, steps (1), (2), and (3) 


A 2 'i+1 “ °* 8 X A 2 + 0 * 2 (x A ' 
"i" denotes the number of the 


i • ^ V Ar ; . 

X 

iteratio 


where the 


ion step. 


When the convergence condition for the point is 

met (usually in two iterations) , the distance x of this 

point becomes the distance for the first iteration in the 

-i+1 -j 

next time step, i.e., x— = (dotted vertical line in 

Fig. 3-3). 


In Fig. 3-4, the lines sent out from A? and A? to the 

i+1 J 

right-boundary point B r (x N ,t J ) are the right-running 
characteristic and particle path, respectively. The calcu- 
lation of the right-boundary values by the simultaneous 
solution of the above-mentioned pertinent system of equations 
is similar to that of the left-boundary values. 

The positions of the points A-^, A 2 , and A^ change 

very slowly with respect to time. The directions of exagger- 
ated shifts shown in Figs. 3-3 and 3-4 are typical for the 
chamber-filling interval of the transient. 


It should be noted that some solutions, based on separate 

19 

integrations along the characteristic lines, can result 
in oscillations of the boundary values. The quasi-linearized 
simultaneous-solution method results in a smooth and stable 
solution. The method may be applied, however, with caution 
to other configurations which cause different boundary con- 
ditions, such as zero aft-end velocity during part of the 
transient due to use of a nozzle closure, and change of 
velocity direction at the fore-end due to short igniter 
operation . 
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3.6 Computation Efficiency and Convergence Tests 

The numerical calculation solves Eq. (3-4) along 
with the analogous equations for T and p, which represent 
the governing equations in their implicit difference algorithm. 
For the N-l spacewise interior points considered, 3(N-1) 
linear algebraic equations have to be solved simultaneously 
for every time step for 3 (N-l) unknowns. The value of N-l 
depends upon the given lengths of the propellant and aft 
sections, x„ - x , and the spatial mesh size. Ax, used 

“ r 

in the calculation. The faster the pressurization process, 
the smaller the value of Ax that must be used, and there- 
fore, the larger the number of interior calculation points. 
Also, the time step, or temporal mesh size. At, should be 
kept small enough, so that the variations in flow properties 
with respect to time can be studied, and the quasi-linearized 
nonlinear and inhomogeneous terms can be properly represented. 

For efficient computation, the finite-difference equa- 
tions are arranged into so-called block-tridiagonal matrix 
24 

form. In this particular form, the 3 (N-l) x 3 (N-l) matrix 

of the coefficients of the unknowns has (3N-5) 3 x 3 square 

sub-matrices as single elements. An economical solution is 

obtained through block-factorization into two block-bidiag- 
. 24 

onal matrices. For accurate computation, the computer 

option of double precision is used in the matrix calculation. 

Variation of the mesh sizes, At and Ax, was carried 

out to test the convergence of the solution. This variation 

was conducted in connection with the so-called Courant- 

Friedrichs-Lewy stability condition for one-step space 

20 

difference equations, which is 

( 1 u | + c) At/ Ax < 1 (3-10) 

However, when an implicit difference scheme is used, the 

20 

numerical solution should be unconditionally stable.* 

Once a stable and uniform solution is obtained, further 
reduction in At does not significantly alter the calculated 

*ln practice, during the initial phases of pressurization, 
the stability conditions for the entrance region and the 
slots require a smaller At than that which satisfies 
Eq. (3-10). 
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results. Mesh sizes, for which the maximum value of (|u| + 
c) At/Ax is slightly greater than one, can often be used. 
However, during the initial portion of the solution much 
smaller values of At are necessary. 

Various tests were performed on the sensitivity of the 
solution to changes in initial conditions and small changes 
in important input parameters, such as the igniter mass flow 
rate, burning rate law, propellant ignition temperature, 
adiabatic flame temperature, and propellant density. In all 
these tests the solutions were bounded and smooth, and changed 
only slightly for small changes in these parameters. This 
demonstrates the existence of neighboring solutions and the 
convergence of the entire numerical solution. 

3.7 Organization of Computer Program 

The program, consisting of over 3000 Fortran 
Statements, is organized into a number of modules (i.e., 
subroutines) which are called from MAIN as indicated in 
Table 3-1. A brief description of each subroutine is given 
in Table 3-2. The program is written in a reasonably general 
form so as to reduce the need for user initiated modifications. 
However, those subroutines which control the input (i.e., 
SDATAC) and the output (i.e., ANS, SRNBUG , THRUSl, PNCH) 
are largely self-explanatory and can be readily interpreted 
by the user. Similarly, those subroutines which calculate 
the interface conditions (i.e., MIGDOT, HCCAL , FCCAL, TSCAL 
and BRCAL) correspond closely to the descriptions given in 
Section 2 and are relatively uncomplicated. Three of the 
subroutines (i.e., TRID, LAGIN and RUNGE) are standard 
packages for performing standard numerical operations. The 
heart of the numerical solution consists of subroutines 
(i.e., PREPP, SOLUTP , PREPC , SOLUTC, INHOMO) coupled to the 
left boundary (LBC) and to the right boundary (RBC) . These 
subroutines should be of little interest to the general user 
and are the domains of persistent numerical analysts. The 
subroutines (i.e., SLOTK, FCT, and FCT2) follow the discussion 
in Section 2 very closely. 



Table 3-1 


FLOW CHART FOR COMPUTER PROGRAM* 



♦Subroutines defined in Table 3-2. 














41 


Table 3-2 

ORGANIZATION OF SUBROUTINES USED IN PROGRAM 

As indicated in Table 3-1, main subroutines are called 
from MAIN, which is the executive portion of the program. 

The functions of the subroutines are indicated by their 
broad classification. 

Initialization of variables and reading of input data: 

SDATA - Initializes variables used in numerical 
solution. 

SDATAC - Reads in data, writes out data, initializes 
physical constants, calculates physical 
properties . 

Computation of parameters at the boundaries of the main 
stream flow: 

MIGDOT - Calculates mass flow rate from igniter. 

TSCAL - Calculates propellant surface temperatures. 

HCCAL - Calculates convective heat transfer coef- 

ficients . 

FCCAL - Calculates friction factors. 

BRCAL - Calculates burning rates. 

LBC - Calculates the left boundary (i.e., igniter 

end) flow parameters. 

RBC - Calculates the right boundary (i.e., nozzle 

end) flow parameters. 

Numerical solution for contributions from slots: 

SLOTK - Solution to p, T, V, and m . for 
each slot. S10t 

FCT2 - Calculates derivatives to be solved by 
RUNGE. Used when pressures in slot and 
main chamber are nearly equal. 

FCT - Calculates derivatives to be solved by 

RUNGE. Used when pressure differential 
between main chamber and slot is significant. 

RUNGE - A fourth order Runge-Kutta solution of the set 
of first order ordinary differential equations 
which described the slot flow. 
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Table 3-2 (Continued) 

ORGANIZATION OF SUBROUTINES USED IN PROGRAM 


Numerical solution of the three PDE's: 


PREPP 


Prepare tables and variables for prediction 
step. 

SOLUTP 

- 

Implicit prediction solution for T, p, 
and u . 

PREPC 

- 

Prepares tables and variables for corrector 
step. 

SOLUTC 

J 

- 

Implicit corrector solution for T, p, 
and u. 

TRID 

— 

Solves the block-tridiagonal matrix of the 
finite-difference equations. Called by 
SOLUTP and SOLUTC . 

INHOMO 


Prepares inhomogeneous terms used in 
implicit solution. 

Time-dependent output: 

ANS 

- 

Executive subroutine that causes the time- 
dependent answers to be written out. 

SRNBUG 

- 

Writes tabular summary of parameters at 
each axial station. 

THRUS1 

- 

Calculates thrust parameters. 

PNCH 

- 

Punches cards that are used as input to 
programs that plot graphs. 

Miscellaneous 

subroutines: 


TIMEST - Sets internal codes whenever time step is 
changed. 

LAGIN - Interpolation used with MIGDOT to obtain 

m ig(t), with RBCAL to obtain r (p) , and with 
SDATA to obtain A (x) and b(x). 

Jr 




4.0 DESCRIPTION OF COMPUTER PROGRAM INPUT 


The discussion in this section coupled with the samples 
described in Sections 5.0 and 6.0 will enable the user to 
prepare the input to the program. 

4.1 Basis for Selection of Dimensional Units 

Since this is a period of transition to the SI 
system of units, and since the computer program employs con- 
ventional propellant configuration dimensions (which in the 
U.S. are specified in inches) and combustion and thermal 
properties (most widely expressed in the CGS system) , an 
effort was made to accommodate the majority of current 
preferences by the following: 

Input: 

. Propellant and motor geometry may be input either 
in centimeters (cm) or inches (in) . 

. Propellant properties are specified in the CGS 

system (modified so that pressure is in atmospheres) 

. Igniter mass flow may be specified in either grams 
per second or pounds per second. 

Output: 

. When propellant and motor geometries are input in 
inches they are converted to centimeters and dis- 
played as part of the initial output. 

. In the time dependent output displays, distance 
along the motor is printed out in both centimeters 
and inches. Pressure is printed out in both 
atmospheres and psia, and thrust is printed out in 
both newtons and pounds force. All other 
outputs are familiar nondimens ional parameters or 
easy to deal with parameters in the CGS system. 

4.2 Definitions of Input Parameters 

The initial input to the program is in four parts: a 

comment card identifying the run, the data by means of NAMELIST, 
table for m^ vs t, table for geometry, and, as an option, 
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burning rate vs pressure. The input and output for typical 
runs are included in Sections 5.0 and 6.0. . 

Each input variable is described on the following pages. 
Unless indicated otherwise by the decimal or lack of decimal 
in the default value, the NAMELIST variables follow the 
usual Fortran conventions on fixed and floating point 
numbers. The usual caution should be used when building up 
a new data set, e.g., 

(1) whenever practical, pattern new data sets after 
sample cases, 

(2) first checkout run should be for a small fraction 
of the ultimate time period, i.e., a short TMAX, 

(3) build up the data decks for complex cases through 
a progression of checkout cases in which the 
individual features are explored separately, 

(4) first few runs should explore the effects of 
increment size on accuracy, e.g., NDELX and DELTAT. 

(5) examine output for instabilities caused by improper 
time and distance increments or by extrapolating 
beyond values in input tables, and 

(6) make sure dependent variable values in input table 
are continuous, i.e., port area versus distance, 
burning perimeter versus distance, and igniter mass 
flux versus time. 

As an option, several of the input parameters can be 
changed during the run by causing NAMELIST NAME to be read 
in at either a prescribed time or nozzle-end stagnation 
pressure. Indeed, causing NAMELIST NAME to be read several 
times during the ignition transient for the purpose of 
increasing At is a necessity if running time is to be re- 
duced . As indicated in the sample problems. At can be 
increased as the ignition and pressurization events progress. 
The user should pay careful attention to learning when and 
how much to increase At. 

thh 


flRTWINAL lb * 
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@ 

Genera l 

Limits, code words, and initial conditions 
Computer Symbol 

Symbol in Text Description 


@uni ts 

CGS ( modified.) 
British 


+ 

TMAX t , Maximum time. When exceeded pro- sec 

max . 

gram reads new data m NAMELIST 

NAME (see page 53) . Should be mul- 
tiple of current value of TPRINT. 

( 0 . 0 ) 


LAMBDA 


TPRINT 

PZONE 


When time increment At is not input, At 
is taken as 1/20 of the hot-gas wave time 
for an increment, Ax. LAMBDA is a factor 
multiplying that At, for the purpose of 
providing a prescribed adjustment to At. 
(1.0) 

Time interval between print out sec 
of answers. 

When stagnation pressure at end atm 
of chamber reaches PZONE, program 
reads new data in NAMELIST NAME. 

(See page 53.) 

(7777.0) 


DELTAT At Time increment, At. When equal to 

zero, program will attempt to cal- 
culate the initial time step. The 
initial time step should be small 
compared to the ratio of the dis- 
tance increment, Ax, to the sonic 
velocity in the igniter gases. As 
the flow field develops. At should 
be rapidly increased using the 
NAMELIST NAME option described on 
page 53. 

( 0 . 0 ) 


In accordance with this period of transition to the SI system 
of units, several of the parameters may be input (as an 
option specified by the input UNIT) in inches and lbm/sec. 
Those inputs are indicated by @ in the right column. All 
other dimensional inputs deal with combustion and thermo- 
dynamics and are input using a modified CGS system. 

■4- t 

Numbers in parentheses are default values built into the 
program. 
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Computer 

Symbol 

NDELX 

TP I 

PAM 

UNIT 


NINERT 

NPNPXT 


@uni ts 

Symbol CGS (modified) 

in Text Description British 

Number of space-wise steps, Ax, 
along the propellant grain. 

(Maximum value is 29) 

(Fixed point number.) 

( 20 ) 

T . Initial temperature in motor K 

P 1 chamber and propellant. Adianatic 

flame temperature is adjusted by 
program whenever TPI ^ TOREF. 

(298.0) 


P 


am 


Ambient absolute pressure. 

( 1 . 0 ) 


atm 


Codeword to indicate dimensional 
units of selected input values. 

If UNIT = -1.0, calories, centi- 
meters, grams, K, seconds and 
atmospheres (i.e., the modified 
CGS system) are used. If UNIT = 
-2.0, those geometry related 
variables indicated by @ in the 
right column are input in inches, 
the igniter mass flow is input in 
lbm/sec, and all other input 
variables are in the CGS system. 
(-1.0) (Not a fixed point number.) 


When NINERT = 1, part of the wetted 
perimeter may be inert (i.e., not 
propellant). When NINERT = 0, all 
the perimeter is propellant. See 
descriptions of geometry input 
table . 

(0) (Fixed point number.) 


When NPNPXT = 1, cards for plot- 
ting p, u, and T are punched 
every TPRINT. 

(0) (Fixed point number.) 
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Motor Configuration Parameters 

The motor geometry input parameters are illustrated 
in Figs. 3-1 and 3-2. The tabular inputs for port area and 
perimeter are described on page 51 . 


@units 


Computer 

Symbol 

Symbol 
in Text 

Description 

CGS (modifi 
British 

AT 

\ 

Throat area. 

cm ^ 
@in 2 

XP 

X P 

Distance downstream at which 
propellant begins. Flow 
region between x = 0 and Xp 
is considered to be inert 
passage with port area Ap . 
Cannot be less than XE/50O* 

TOT 

cm 

@in 

XG 

X G 

Position at aft end of propellant 
grain. Set by program as 
x G = x E - Ax/2. 

cm 

,@in 

XE 

X E 

Position at end of flow passage 
(XE > XG) and before isentropic 
nozzle flow begins. Distance 
increment is (x E - Xp)/NDELX. 

cm 

§in 


Combustion 

Gas Properties 


GAMA 

Y 

Effective ratio of specific 
heats of all combustion gases. 


W 

W 

Average molecular weight of 
all combustion gases. 

g/g-mole 

TIGN 

T. 

ig 

Effective mean temperature of 
the igniter gases. 

K 

TFREF 

T f,ref 

Adiabatic flame temperature at 
reference conditions, i.e., 

Pref / T 0,ref • (Actual value 
should be used rather than the- 
oretical value. TFREF should be 
consistent with input value of 
W and GAMA and measured c* values 
See discussion in Appendix A.) 

K 

• 
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Friction 

Factor and Heat Transfer Coefficient 


Computer 

Symbol 

Symbol 
in Text 

Description 

@units 

CGS (modified) 
British 

RUFSUR 

e s 

Roughness of the port walls 
(used prior to ignition) . 
When zero, friction is zero. 

cm 


( 0 . 001 ) 

w 

DDRG A factor multiplying the friction 

■ factor, for the purpose of providing 
a prescribed adjustment. 

( 1 . 0 ) 

DDHC A factor multiplying the heat transfer 

coefficient, for the purpose of provid- 
ing a prescribed adjustment. 

( 1 . 0 ) 


Propellant Properties 


FKPR 

X pr 

Thermal conductivity of solid 
propellant. (Btu/sec-in-°R = 
178.6 cal/sec-cm-K) 

cal/sec-cm-K 

ROPR 

p 

K pr 

Density of solid propellant. 
(1 lbm/in 3 = 27.68 g/cm 3 ) 

g/cm 3 

CPR 

C pr 

Specific heat of solid propellant. 

cal/g-K 

TOREF 

T 0,ref 

Initial propellant temperature 
for reference conditions. 
(298.0) 

K 

SIGP 

° P 

Temperature sensitivity of K 

burning rate at constant pressure, 

(3 In r/3T n ) . Used only with Eq. (4-1). 
(0.002) u p 

TPSCRI 

T ps,ig 

Surface temperature at which 
propellant ignition occurs. 
(600.0) 

K 


Burning 

Rate Relationship * 



The nonerosive burning rate, r q, can be specified 
either by the empirical equation 

r 0 = r r e f <p/£> ref )n 8 x P Io p (T pi ' T 0,re£ n (4 ' 1) 

*¥fMF,""and"RREP form a consistent set of parameters evalu- 
ated at a prescribed reference condition, (e.g., »pOREF = 
298 . OK, PREF = 68.08 atm). 
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or by a table (which does not include c correction) : 

, P 

r Q = f (p) (4-2) 


The following inputs apply to the empirical equation and the 
erosive burning contributions. See page 52 for the inputs 


that prescribe the 

r = f (p) table. 


Computer 

Symbol 

Symbol 
in Text 

Description 

@units 

CGS (modified) 
British 

RREF 

r ref 

Burning rate at reference con- 
ditions. If NDATA > 0, steady 
state burning rates are obtained 
from RDATA vs PDATA table. 

cm/sec 

PREF 

^ref 

Absolute pressure for reference 
condition propellant properties. 
(68.08) 

atm 

BREXP 

n 

Exponent in burning rate 
relationship. 



If the following modified Lenoir-Robillard erosive burning 
rate law is used: 


r = r n + k h exp(-8 rp /G) (4-3) 

0 e c e pr 

where G is (pu) , and h is the convective heat transfer 
coefficient, then c 

3 

EBC k Pre-exponential constant in cm -K/cal 

' Lenoir-Robillard erosive burning 

relationship, Eq. (4-3). 

(C.O) 

EBEX 8 Factor in exponent of Lenoir- 

8 Robillard erosive burning 

relationship, Eq. (4-3). 

(125.0) 

If the more conventional form of the Lenoir-Robillard erosive 
burning rate law is used 

r = r 0 + a e G °* 8d h~°* 2 ex P (- i e e r P pJ /Gl) (4-4) 

where hydraulic diameter is d^ = 4A /b, then 
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EBC a 

e 

EBEX -$ 

e 

To convert a 

e 

by 0.102. 


Pre-exponential factor in Eq. (4-4) . 

Code word and factor in exponent 
of Eq. (4-4). When < 0, signal 
to use Eq. (4-4) . 

from in^ * ®/lbm^ ' ®-sec^‘ ^ to CGS, multiply 


Parameters- Used in Thrust Calculation 


DE 


CM 


ALFAD 


For details of thrust calculation, see Appendix 1. 

D Exit plane diameter. If 0.0, cm 

thrust calculation is bypassed. @in 

( 0 . 0 ) 

c Nozzle coefficient for thrust 

losses. 

( 0 . 0 ) 

Nozzle divergence 1/2 angle. deg 

N (15.0 


EROAT E t 

EROEXP n t 


Parameters in empirical nozzle cm/sec 
erosion equation, dD/dt = 2* EROAT @in/sec 
P**EROEXP . 

( 0.0 & 0 . 8 ) 


Igniter Mass discharge Versus Time Table 

First table following NAMELIST NAME input. 

NIGTAB* Number of points in table. (Max- 

imum is 30) 

Fixed point number; part of 
NAMELIST NAME. 


Use 2E10.0 Format for TIGTAB and MIGTAB. 


Two point interpolation is used. 


TIGTAB 

t 

Time after beginning of run 
Second time in table should 
at least 30 At's. 

be 

sec 

MIGTAB 

m. 

ig 

Mass flow rate from igniter 
Program will assign a small 
to m^g at t = 0 so as to 
isfy starting conditions of 
solution. 

• 

value 

sat- 

g/sec 

@lbm/sec 


•include as part of NAMELIST NAME input. 
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Propellant Geometry Input Table 

Second table following NAMELIST NAME input. For 
more details see Figs. 3-1 and 3-2. 

Input which apply to both monolithic and segmented motors: 


Computer 

Symbol 

Symbol 
in Text 

description 

@units 

CGS (modified) 
British. 

NAPDVX* 


Number of points in table and 
stations.* (Maximum is 30) 
Fixed point number, part of 
NAMELIST NAME. 



Use 8 El 0.0 format for XAP, APD, BPD, WPD, XSLOT, DELABS, 

DELVS and WIDTH. Values in table must be consistent with 
XP, XG and XE. 

When circumferential slots are considered, the table 
must include one set of entries for each x station corres- 
ponding to equally spaced stations Ax apart. The x values 
are to be specified so as to correspond to the leading 
edge of each increment. 


XAP 

X 

Distance from head end of motor. 

cm 

@in 

APD 

A P 

Port area of main channel. 

_ 2 
cm 

@in 

BPD 

b 

Perimeter of propellant surface 
in main channel . 

cm 

@in 

WPD 

Pi 

Total perimeter of propellant and 
inert surface (i.e., wetted 
perimeter) before ignition. 

cm 

@in 

Inputs which 

apply only to segmented motors:* 


XSLOT 

X slot 

Distance to the upstream tip of 
the slot from the head end of 
the motor. 

cm 

@in 

DELABS 

5A b,slof 

Burning surface area in slot. 
(Subtract bw s i Q t from actual 
burning surface area in the slot.) 

2 

cm 

@in 


♦When segmented motors are considered, input geometry for each 


station, i.e., NAPDVX = NDELX + 1. 
inputs corresponding to slot at x 


When x 
X slot 


< x. 


lot 


< x + Ax, 


are prescribed, 



Computer' 

Symbol 

DELVS 


WIDTH 


Symbol 
in Text 


<$V 


slot 


w slot 


•52- 


aEP KODUCBILnY O^HE 

■sfiimNAL PAGE IS W 


Description 


&unit8 

CGS (modified) 
British 


, 3 

Volume (region of gas not in the cm 3 

flow field) of slot. Program @in J 

treats this as volume sink and 
not part of the flow passage. 


Width of axial separation cm 

between segments. (Must be @in 

sufficiently large so that 
flow out of slot is less than 
about Mach =0.3)* 


Burning Bate Table (r vs p) 


Optional third table following NAMELIST NAME. Not 
used when r = r re f (P/P re f ) n is desired. 

Use 2E1Q.0 format for PDATA (in columns 1 to 10) and RDATA 
(in columns 11 to 20) . Table is logged by the program, and 
interpolation is in In r vs In p table. 

NDATA* * Length of table (i.e., the number 

of r,p pairs) . Maximum table 
length is 30. If zero, table not 
used. 

( 0 ) 


NPPR** Number of points to be used in 

interpolation. Linear interpolation 
is prescribed by inputing 2. If 
three or more points are used it is 
important that the values in the tables 
be uniform and continuous. 

( 2 ) 

PDATA Values of pressure in r vs p atm 

table. (Place in columns 1 to 10) 

RDATA Values of burning rate, r cor- cm/sec 

responding to pressure, PDATA. 

A point in the table must corres- 
pond to PREF and RREF. (Place in 
columns 11 to 20) 


♦Maintaining low Mach numbers out of the slot is a practical 
design consideration, which is discussed in Ref. 27. 

**Read in as part of NAMELIST NAME. 



4.3 Intermediate Insertion of New Input Parameters 

Ther'e are various situations where it is either 
convenient or necessary to interrupt the calculations and 
insert new data. For example, to decrease the running time, 
the time step. At, must be increased as the flow develops. 

To accomplish this, the program will read in the NAMELIST 
NAME when nozzle-end stagnation pressure > PZONE or 
time > TMAX. This is illustrated in Tables 5.1 and 6.2. 

The input parameters which can be read in include: 

Symbol 

DELFAC (new value of At = DELFAC. *At) 

TMAX (should be multiple of current value of TPRINT) 
DELTAT (new time step, At) 

TPRINT 

AT (throat area) 

PZONE* NOTE * 

TIGN 
DDRG 
DDHC 

For example, if at t = 0.01 sec the following changes are 
desired: 

(1) increase At by 40%, 

(2) print out answers every 0.005 sec, 

(3) next call to read NAMELIST NAME is to be at 
either t = 0.04 sec or p E gtag = 2.0 atm, 

the input should be arranged as follows: 

(1) original input set should contain TMAX * 0.01. 

(2) after geometry table place input array for NAME- 
LIST NAME , &NAME DELFAC = 1.4, TPRINT = 0.005, 
TMAX = 0.04, PZONE =2.0 SEND 

♦Because of internal pressure units, this value must be input 
in atm and not in psi. 


Unite 


sec 

sec 

sec 

2 

cm 

atm 

K 
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4,4 Punched Output for Input to Plot Programs 

As specified by the codeword NPNPXT, punched cards 
containing many of the output variables can be obtained. 
Using this punched output, the user can prepare a special 
purpose program to plot the results. 

At the beginning of each punched output deck is the 
title card that goes with the data set. The following 
values are punched in the next card using a 12A4 field: 



Distance downstream at which cm 

propellant begins. 


x_ Position at end of flow passage cm 

and before isentropic nozzle 
flow begins. 

2 

A fc Throat area. cm 

BLANK 

Initial gas temperature in K 

chamber . 



Pint 

M int 

Ax 


Initial chamber pressure. atm 

Initial Mach number 

Distance between axial stations. cm 


UNIT Codeword to indicate units used sec 

in program. 

t Time punched data begins. sec 

y Ratio of specific heats. 

R J/W Gas constant divided by average 

u molecular weight. 

The following sets of values are punched for each TPRINT 
interval. The first card in the set contains (in an E14.6, 
218, 2F10.1 field) : 

t Time sec 


NUT 


Total Number of time steps used 
to this point. 
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NCARDS Number of cards corresponding to 

each time interval. 

P E stag Nozzle end stagnation pressure. atm 

F Thrust newtons 

The next NCARDS cards in the set contain NDELX+1 pressures, 
followed by NDELX+1 temperatures, and followed by NDELX+1 
Mach numbers (using a 20A4field) . 
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5.0 SOLID ROCKET BOOSTER CONSIDERED AS MONOLITHIC CONFIGURATION 

The results presented in this section along with those 
of the next section serve two purposes: (1) they define 

and demonstrate the use of the computer program and (2) 
they provide information which is applicable to SRB’s. The 
calculations presented in this section were performed during 
the time the procedures for analyzing the slot interactions 
were being developed. Accordingly, in this section the SRB 
is treated as a monolithic motor with a total chamber volume, 
length, and burning surface area similar to that of the 
segmented motor described in Ref. 25. Since the scope of 
the study did not permit all of the parametric studies to be 
redone using the final version of the computer program and 
since the results are generally applicable to the SRB 
described in Ref. 25, they are included in this report. In 
Section 6, the final version of the computer program is used 
to analyze the SRB as a segmented motor. 

The main approximation in using a monolithic motor 

configuration to analyze a segmented motor is that the 

segmented motor geometry must be approximated in terms of 

port areas and perimeters. However, since the burning 

surface area and chamber volume in the slots of the SRB 

25 

under development are relatively small, the approximation 
is reasonable. 

The input prepared to simulate the SRB design is shown 
in Tables 5-1 and 5-2. Table 5-1 is a listing of the input 
data cards and Table 5-2 is a copy of the printout prepared 
by the program to display the input. Note that the variables 
in Tables 5-1 and 5-2 are the same as described in Section 4. 
Preparing the input for the complex head-end geometry of the 
Ref. 25 SRB requires careful attention. It is important to 
account for all of the burning surface area, chamber volume, 
and port area. All of this must be done using a few dis- 
tance increments and without allowing an abrupt discontin- 
uity in the flow field. The reader should pay particular 
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Table 5-1 

INPUT DATA FOR MONOLITHIC MOTOR (SRBl) 

S9B1 TPSCRIR50 KOEPO 24 

SHAHE 

T*AX=.002 
r "?PIN’ 1, = .00 2 
DELTA ,r = .00016 
HDETX=21 


0HlT=-2. 
HPSPXT=0 
AT=2327. 
XP-3 . 

XR® 134 1.3 
XE=1377. 

GAH A=1 . 136 3 
8=28.21 
TFBEE=3361. 
TTGN=?650. 

TtnFsn”=o. oi 


NOTE : LEAVE COLUMN 

1 BLANK IN ALL 
NAMELIST INPUT 


FFPP=0. 00 1 1 
POPP = 1.758 
C?R=.3 
TPSCST=850. 
RPEE=1 .0783 
BEEXP = 0.35 
EBC=0. 
BE=145.64 


CH= .98 
ALFAD=12. 31 
NTG' T ’AB=6 


vap7)vx=io 



5 EH D 




0. 

25. 



.020 

560. 



.250 

560. 



.340 

300. 



.450 

100. 



.50 

0. 



0. 

2554. 

1018.92 

1018.92 

3. 

2554. 

1018.92 

1018.92 

68. 4286 

2554. 

1018.92 

ini8.92 

133.8571 

255 4. 

1818.92 

1018.92 

199. 2857 

2870. 

611.23 

611. 23 

264.7143 

3187. 

203.53 

203. 53 

941.6 

3187. 

203.53 

203.53 

1023. 

3096. 

203.53 

203.53 

1286. 

3400. 

, 203.53 

' 203.53 

1377. 

4500. 

203.53 

203.53 

6HAHE 




TMAX= 

.01 




81L?AC=1. 3 


KERB 


TMAX=. 
THAX=, 
DEIFAC=2. 


151 
12 J 


Last value is 
retained in 
memory 



TH AX=.2 
<!>PR-|J'T’=. 005 


TH AX= . 4 
TPFINT=.01 


SEND 

6HAME 


TH A X=. 5 
THAX=.6 



REPRODUCIBILITY OF TEE 
"ORKKNAL PAGE IS POOP 
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Table 5-2 MONOLITHIC INPUT TO SIMULATE SRB DESIGN 

HIGH PERFORMANCE SOCKET MOTOR IGNITION TRANSIENT PREDICTION PROG RAH 
INCLUDES : SPATIAL S TIME DEVELOPMENT OB’ P, 0, ST AND FLAHF SPREADING 


SRBl 24 

GENERAL: LIMITS, CODE WORDS, £ INITIAL CONDITIONS: 


TH AX 

LAMBDA 

TPSINT 

PZCNE 

DELTAT 

NDELX 

: 0.002 

1 .00 

0.001 

7777. 

00 0.0001600 

21 

TPI 

PAM 

UNIT 

NINER 

T NPNPXT 


298.0 

1.C0 

-2 . 

0 

0 


MOTOR CONFIGURATION PARAMETERS: 



AT 

XP 

XG 

XE 



2327.00 

3.00 

1374.00 

1374.00 



COMBUSTION GAS 

PROPERTIES 




GAMA 

w 

TERR? 

TIGN 



1. 136 

28.21 

3361.0 

2650.0 



FRICTION 

TRANS? 

ER AND HEAT 

TRANSFER 

COEF. : 


RUFSUR 

DDRG 

. DUHC 




0. 0100 

1.00 

1 . 00 




PROPELLANT PROPERTIES: 




FKPR 

ROP* 

CPR 

TOREF 

SIGP 

TPSCRI 

0.00110 

1.758 

0. 30 

298.00 

0.CC200 

850.00 

BURNING 

RATE RELATIONSHIP: 




FREE 

PR EB 

UREXP 

EBC 

FEEX 


1 .078 

6P.08 

0. 3500 

0.0 

125.000 



PARAMETERS USED TN 

THRUST 

CALCUI ATICNS : 

DE 

CM 

ALB- AD 

EFCAT ERCIXP 

145.64 

0.980 

12. 11 

0.0 O.P0 

IGNITER 

HASS DISCHARGE VS. 

TIME TABLE: 

TIME 

FLOW PATE 



0.0 

25.0 



0.0200 

560.0 



C .2500 

560.0 



C. 3400 

300.0 



0.4500 

100.0 



C.5000 

0.0 




PROPELLANT GEOMETRY TABLE: 


LENGTH 

PORT AREA 

PROPELLANT 

WETTED 



PERIMETER 

PERIMETER 

0.0 

2554.00 

1018.92 

1018.92 

3.00 

2554.00 

1018.92 

1018,92 

68.43 

2554.00 

1018.92 

1018.92 

133. 86 

2554.00 

1018.92 

1018.92 

199.29 

2870.00 

61 1.23 

611.23 

264.71 

3187.00 

203.53 

203.53 

961.60 

3187.00 

203.53 

203.53 

1023.00 

3096.00 

203.53 

203.53 

1286. OC 

5595. 00 

203.53 

203.53 

1377. 00 

5595.00 

203.53 

203.53 
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attention to entries at stations 4, 5, and 6 in Table 5-1 and 
note that station 5 is a transition section. 

A typical time dependent output is shown on Table 5-3. 

The output has been annotated to define the output values. 
Figure 5-1 shows the calculated pressure-versus-time at four 
positions along the port. Also shown is an expanded portion 
of the curve emphasizing the reflection of the first strong 
wave (that arrives at the nozzle closure) and the times at 
which ignition occurs. Immediately apparent are the inter- 
actions produced by the nonuniform pressure at 0.08 seconds 
which are a result of the high velocity gases being driven 
ahead of the ignition front and reflecting off the nozzle 
converging section. (It is important to note that this 
nonuniform pressure is not expected to produce any structural 
problems since it occurs at an intermediate pressure.) The 
nonuniform pressure is to some degree a characteristic of 
this particular type of motor, i.e., a long motor with a 
large burning surface area in the head end and a moderately 
small port- to- throat area ratio. Figure 5-2 shows some of 
the details of the, pressure wave development in the p-x 
plane and clearly indicates the nature of the wave reflection. 
When the results were tested by decreasing the time and 
distance steps, no important differences in the wave character 
or magnitude were noted. 

The results of Fig. 5-1 should not be considered as 
detailed predictions of the SRB since much is yet to be 
learned from correlating previous large motor firing data. 

In particular, users will have to gain experience with pre- 
dicting convective heat transfer rates and propellant igni- 
tion criteria. 

The ripples on the p vs t curves on Fig. 5-1 are a 
result of longitudinal waves reflecting off the nozzle and 
head end closures. Note that these oscillations damp as the 
motor reaches its operating pressure. The oscillations are 
a physical characteristic of this particular SRB (as opposed 
to an indication of numerical difficulties) and are not 
observed in all calculations of large booster ignition 
transients . 



Table 5-3 (Part 1) 

OUTPUT FROM MONOLITHIC MOTOR RUN (SRB1) 
- DURING FLAME SPREADING - 


TIME -0.8400E-01 ( 196 DEL TAT » DEL TAT =0.5000E-03 


I* DISTANCE PRESSURE T MACK RATE BfiTORN TAUB FRICT. CftNV. M>A TSURF PORT 



CM 

IN 

AIM 

PS I A 

DEO K 

NO. 

CM/ SEC 

CM FACTOR COEF • G/SEC CMSO 

K 

CMSO 

6 

0.0 ' 

0.0 

14.52 

213.3 

2654.8 

0.086 

6.63 

1.000 0.01 TO 

0.0 

0.0313 

15.09 

850, 

16521*4 

* 

173.8 

68.4 

13.86 

202.8 

2969.7 

0.240 

0.62 

1.000 6.0123 6.0 

0*0613 

38.39 

850. 

16509.2 

M 

340.0 

133.9 

12. 15 

178.6 

3103.3 

6.426 

0.69- - 

1.000.0.0094 0.0 

0.0827 

58,32 850* 

16601*7 

3 

506.2 

199.3 

9.46 

139.0 

3114.5 

0.644 

0.54 

1.000 0.0067 

0.0 

0.0859 

68.70 

850. 

18526.4 

4 

672.4 

264.7 

7.81 

114.8 

3096. 1 

0.788 

0.51 

1.000 0.0045 

0.0 

0.0750 

69.64 

850. 

20563.5 

6 

838.6 

330.1 

7.55 

110.9 

3143.4 

0.802 

0.50 

1.000 0.0036 

0.0 

0.0723 

67.95 

850. 

20563.0 

b 

1004 .8 

395.6 

7.30 

107.2 

3183.3 

0.812 

0.49 

1.000 0.0028 0.0 

0.0698 

66.06 

850. 

20562.7 

r 

1116.9 

461.0 

7.01 

103.6 

3190.3 

0.833 

0.49 

i.eoo 0.0021 

0*0 

0.0680 

- 65*05 

850* 

26562.3 

is 

1337.1 

526.4 

6.73 

98.9 

3130.6 

0.873 

6.48 . 

1*600 0.0012 0.0 

0.0678 

66,06 

850, 

20561.8 

9 

1503.3 

591.9 

6.41 

94.1 

£983.7 

0.916 

0.47 

1.000 0.0005 

0.0 

0.0675 

67.57 

850. 

20561 .5 

10 

1669.5 

657.3 

6. 18 

90.7 

2807.6 

0.931 

0.0 

0.0 0.0 

0.0031 

0.0667 

68.29 

828. 

20561.2 

11 

1835.7 

722.7 

6.07 

89.2 

2654.8 

0.922 

0.0 

0.0 0.0 

0.0031 

0.0657 

68.33 

775. 

20561.2 

12 

2001.9 

788. 1 

S.97 

87.7 

2514.7 

0.909 

0.0 

0.0 0.0 

O .0030 

0.0647 

66 4 06 

721. 

20561 .2 

13 

2168.1 

653.6 

5.88 

86.3 

2363.5 

0.898 

0.0 

0*0-- 040: 

040030 

0.064 1 

68,34 

663. 

20561.2 

t* 

8334*3 

919*0 

5.77 

84.8 

2164.6 

0.896 

0.0 _ 

0.0 6.0 

0.0030 

0«064^ 

69,94 

601. 

20561.2 

lb 

2500.4 

984.4 

5.55 

81.6 

1895.4 

0.906 

0.0 

0.0 0.0 

0.0030 

0.0649 

72.82 

536. 

20342.9 

16 

2666.6 

1049.9 

5. 11 

75.0 

1569.4 

0.923 

0.0 

0.0 0.0 

0.0029 

0.0640 

75.01 

472. 

20174.4 

IT 

2832.8 

1115.3 

4.54 

66.6 

1230.8 

0.919 

0.0 

0.0 O.G 

0 . 0029 

0.0605 

75.13 

413. 

20662.3 

16 

Z999.C 

1 180.7 

4.02 

59.1 

929.9 

0.879 

0.0 

-■ .0*0:'/ ■ . 

Q.0O29 

0.0660 

73.40 

369. 

21150.3 

» 

3X65.2 

1246.1 

3.66 

53.8 

694.9 

0.780 

0.0 

6.0- 0*6. 

0.0029 

0.0476 

68466 

339, 

21638.2 

as 3331.4 

1311.6 

3.65 

53.6 

526.0 

0.574 

0*0 

0.0 0.0 

6.0029 

0.0373 

68*11 

318. 

23929.6 

21 

3497.6 

1377.0 

4.22 

61.9 

408.1 

0.326 

0.0 

0.0 0.0 

0.0030 

0.0 

42.86 

306. 

29032.2 

PESTAG.ATM= 4 

.48 PESTAfa.PSIA= < 

85.79 

MASS BURNED- 0.28148 I 

05 






THRUST 

THRUST CELM 

AT 


EPSGNG 

EPSON PA/PE 






NEWTONS 

LBF 



CH**2 


— 

- ■ 






663260.3 

149112.0 

0.9738 15012 

.867 

7.1590 

2.0637 

1.783 






Table 5-3 (Part 2)* 


OUTPUT FROM MONOLITHIC MOTOR RUN (SRBl) 
- DURING RAPID PRESSURE RISE - 


TIME -0.13006 O0 t 2*8 KLUT) 


n 


X 

P 



M 

r 

r/r 0 f 

N 

DISTANCE 

PRESSURE 

T 

MACH 

RATE 

BR/bRN TAUB FR1CT. 


CP 

IN 

ATM 

PS1A 

DEO K 

Nli. 

CM/SEC 

CM FACTOR 

o 

0.0 

0.0 

23«6? 

347.7 

2653.1 

0,062 

0.74 


1 

173.8 

68.4 

22.96 

337 .3 

2944.1 

0.176 

0.74 

1.006 0.0427 6.0 

2: 

340.0 

133.9 

21.47 

315.4 

3077.3 

0.315 

0.72 

1.00 0 0.0 398 0.0 

3 

506.2 

199.3 

19.08 

280.3 

3105.8 

0.435 

0.69 

1.000 0.0353 0.0 

4 

672.4 

264.7 

18.77 

275.7 

3122.8 

0.427 

0.69 

1.000 0,0322 0.0 

6 

839.6 

330.1 

19.35 

284.3 

3156.0 

0.400 

0.69 

1*000 0.0314 0.0 

6 

1004 .9 

395.6 

19.20 

282.0 

3176.6 

0*409 

0.69 

1.000 0.0306 0.0 

T 

1178.9 

461.0 

19.65 

288.7 

32X1.3 

0.358 

0.70 

1.000 0.0297 0.0 

8 

1337.1 

526.4 

20.26 

297.6 

3249.8 

0.360 

0.71 

1.000 0.0286 0.0 

9 

1503.3 

591.9 

21.19 

31 1 .3 

3294.4 

0.320 

0.72 

1.000 0.0278 0.0 

10 

1669.5 

657.3 

22.07 

324 .3 

3338.0 

0.283 

0.73 

1.000 0.0271 0.0 

11 

1835.7 

722.7 

22.52 

330.9 

3373.3 

0.267 

0.73 

1.000 0.0264 0.0 . 

12 

2001.9 

788.1 

22.35 

328 .5 

3396.9 

0.273 

O. 73 

1.000 0.0257 0.0 

13 

2168.1 

853^ 

21*92 

322*0 

3414.2 

0.295 

0.73 

1.000 0.0250 6*0 

14 

2334.3 

919.0 

21.35 

313.7 

3428.2 

0.316 

0.72 

1.000 0.0240 0.0 

15 

2500.4 

984.4 

20.93 

307.4 

3443.5 

0.343 

0.71 

1 .OOO 0.0233 0.0 

16 

if666 «»6 

1049.9 

20.30 

298,2 

3454.2 

0.365 

0.71 

1.000 0.0226 0.0 

17 

2832.8 

1115.3 

19.94 

292*9 

3469.7 

0.366 

^0*70 

• 1.660 0.0217 0.0 r 

la 

2999.0 

1180.7 

19.38 

284.8 

34B1.3 

6^91 

0*69 

. 1* 600 - 0^02l 4^ ^ 0it^ ~ -rrT 

19 

3165.2 

1246*1 

19.22 

282*3 

3499.3 0.411 

6-69 

1.060- 0.0202 -0.0 

20 

3331.4 

1311.6 

18.76 

275.5 

3514.6 

0.371 

0.69 

1.000 0.0154 0.0 

21 

3497.6 

1377.0 

19.52 

286.8 

3524.2 

0.326 

0.0 

0.0 0.0 0.0 


At 

OlOLTAT *®»3 


TSURF 
3 K 


PORT 

CMSQ 


0.6314 

o.ono 

0.1027 

0.1099 

0.0925 

0.0881 

0.0075 

0.0042 

0.0801 

0.0746 

0.0691 

0.0663 

0.0665 

0.0680 

0.0708 

0.0740 

0.0764 

6.9770 

0.6754 

- ■ V-44-4 4 w 9 1 . 

0.0688 

0.0 


. 15.13650. i 6603. 3 
: 46.04 850. l*59ftiT 
76.65 850. 16580.2 
93.61 850. 18570.7 
90.25 850. 20577.7 
66.76 850. 2657T.3. , 
67.59 850. 26*76.6 o> 


80.53 850. 20575.9, 
74.32 650. 20575.4 
68.10. 850. 20575.1 
66..T* 656. 20674.7 

- - 4MLU i 

- 97174 - K l|g-f 6B -6-. 

60.56 8SO. 29674.0 
72.47 850. 20573.5 
77.01 850. 20354.9 
79.35 650. 20186.0 


73.79 850. 23937.5 
67.51 660. 29040.0 



CfXm 

<7i_m 


A t 

AT . . 

CM»*2 



945365.8 


1.3335 15012.867 


ex 

7.1590 


ex 

7. 1590 


2.105 


>ut symbols identified by annotations. 


Table 5-3 (Part 3) 


OUTPUT FROM MONOLITHIC MOTOR RUN (SRBl) 
- APPROACHING FULL CHAMBER PRESSURE - 


L . — ‘ — — 



■■ - • 



TIME =0.4200E OO 

i 544 OELTAT1 DELTA! = 

T 

N DISTANCE 

PRESSURE 

T 

MACH 

RATE 

BR/BRN 

TAUB ERICT. 

CONV. 

M/A TSURF 


JIN. • - 

ATM 

PSlA 

DEG It 

NO. 

CM/SEt 


CM FACTOR 

COFF. g/SEC CM5G 

X 

O 0.0 

0.0 

59.43 

873.0 

2651.1 

0.006 

1.03 

1.000 

0.3271 

0.0 

0.0109 

4.03 

850. 

1 173.8 

88.4 

58.90 

865.3 

2924.6 

0.076 

1 .03 

1 .COO 

0.3208 

0.0 

0.0780 

52.33 

850. 

2 340.8 

133.9 

58.07 

853.1 

3055.8 

O. 152 

1.02 

1 . GOO 

0.3149 

0.0 

0.1265 

100.22 

850. 

- - *»*■ -- M w* #jgL- . .. 

s|S?S3 

56. 08 

823.8 

3097.3 

0.207 

1.01 

1.000 

0 .3059 

0.0 

0.1434 

131.06 

850. 


284.7 

55.78 

819.4 

3110.7 

0.208 

l.ol 

1.000 

043023 

0.0 

0.1242 

130.84 

850. 

"•S' 838.6 

330.1 

55.71 

818.4 

3123.4 

0.217 

1.01 

1.000 

0.3015 

0.0 

0.1256 

135.84 

850. 

6 1004.8 

396.6 

55.21 

81 1 -O 

3132.5 

0.232 

1.00 

1.000 

0.2997 

0.0 

0.1293 

143.71 

850. 

7 1170.9 

461.0 

54.95 

807.3 

3142.9 

0.244 

1 .00 

1.000 

0.2985 

0.0 

0.1324 

150.39 

850. 

8 1337.1 

626.4 

54.48 

800.2 

3151.2 

0.258 

1.00 

1.000 

0.2966 

0.0 

0.1356 

1 57 . 27 

850. 

9 1503.3 

591.9 

54. 16 

795.4 

3159.8 

0.272 

1.00 

14000 

0.2952 

o.© ; 

O. 1392 

164.63 

850. 

TO 1669.5 

657.3 


TB7.7 

3166.5 

0.285 

0.99 

i .000 

0.2935 

0.0 

0.1420 

170.85 


111835. 7 

722.7 

53.26 

782.4 

3173.7 

0.301 

0.99 

1.000 

0.2920 

0.0 

0.1461 

178.80 

.850. 

12 2001.9 

788.1 

52.64 

773.3 

3176.6 

0.314 

0.99 

1.000 

0.2901 

0.0 

0.1484 

184.20 

850. 

13 2168.1 

863.6 

52.29 

768.1 

3184.8 

0.331 

0.98 

i.ooo 

0.2886 

0.0 

0.1529 

192.96 

850. 

14 2334.3 

919.0 

51.53 

767.0 

3187.8 

0.343 

0.98 

1 .GOO 

0.2863 

0.0 

0.1545 

197.17 

850. 

1*2*00.4 

984.4 

51.19 

752.0 

3193.3 

0.367 

0.98 

1.000 

0.2847 

0.0 

0.1612 

209.23 

850. 

I&2685.6 

1049.9 

50*10 

736.0 

3193.2 

0.382 

0.97 

1.000 

0*^S2O 

0.0 

0.1626 

212.97 

850. 

17 2832. B 

1115.3 

50.01 

734.8 

319928 

0.398 

0.97 

1 .009 

0.2807 

0.0 

0.1664 

221.20 

850. 

18 2999.0 

1180.7 

49. 08 

721 . 1 

3200. 7 

0.393 

0.96 

1 . ooo 

0.2785 

0.0 

0. 1612 

214.72 

850. 

19 3165.2 

1246.1 

49.39 

725.5 

3209.2 

0.411 

0.96 

1.000 

0.2780 

0.0 

0.1665 

225.59 

850. 

20 3331.4 

1311 .6 

48.58 

713.7 

32 1 1 . 8 

0.364 

0.96 

1.000 

0.2717 

0.0 

0.1468 196.25 

850. 

UBI 1377.0 

50.99 

;i 7 ^V2r-323 1*4 

04325 

0.0 

~y 

~0*O 

040 

. 040 -rfg 

:\-*a*40* 

850. 


CMS6 


. 17323.1 
. 17306.8 

. 17291.6 

. 18990.2 
. 20716.8 
. 20716.4 

• 20716.5 

. 20714.9 ^ 

. 20713.9 ^ 
. 20713.1 
. 20712.3 
. 20711 .5 ; 

. 20710.6 
. 20709*833 
. 2070£gfcW 
» 2<)489^is^j 
. 2031^60 
. 20806% ; ! 

. 2129^^ 

. 2178^3^ 

• 2406^^1 
. 29172* ti 


54.03 PEST16,P5U* 793.69 MASS BURNED* 0.1 50 7£ 07 


THRUST THRUST 
NEKTONS LBF 


CFtN 


AT 

CM**2 


EPSONG EPSON PA/PE 


U5«M 


PRESSURE, ATMOSPHERES 

OO 10-00 20.00 30-00 4:0.00 50.00 60.00 70.00 
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IGNITION 



FIG. 5-1 PRESSURE DEVELOPMENT AT FOUR STATIONS ALONG THE PORT SHOWING 

THE REFLECTION OF THE GASES DRIVEN AHEAD OF THE IGNITION FRONT- 
(SRB1) 



PRESSURE, fiTMOSPHERES (X10 2 

.00 0.10 0.20 0.30 0.140 0.50 0.60 
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V IGNITION TIME 



FIG, 5-2 PRESSURE WAVE PROPAGATING DOWN THE PORT AND REFLECTING OFF 
THE AFT CLOSURE 
(SRB1) 
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As shown in Figs. 5-3, 5-4, and 5-5, the developing 
flow field during ignition is influenced greatly by the heat 
loss to the propellant. As the gases are accelerated their 
static temperature decreases and as they heat the surface 
their stagnation temperature decreases, which in turn 
decreases their linear velocity. In Figs. 5-3, 5-4, and 5-5, 
the first curve (at t =0.04 sec) is before first ignition 
has occurred; the next three curves are during flame spread- 
ing (the point of flame spreading is indicated by a triangle) . 
As the pressure increases following full ignition, the flow 
field becomes more fully developed, gas velocity decreases 
as the gas becomes more dense, and the gas temperature 
becomes more uniform. The last curve (at t = 0.5 sec) 
approximates fully developed flow conditions and is beginning 
to take on the quasi-steady flow characteristics which are 
adequately represented by p(x) models. 

Flame spreading rates for a range of input conditions 
are shown on Fig. 5-6. The heavy dashed curve corresponds 
to SRBl whose distinguishing characteristic is that the 
ignition temperature is 850K rather than 750K used in most 
of the calculations. 

After the first upstream portion of the grain is ignited, 
the mass flow generated by the main propellant grain soon 
exceeds the mass flow of the igniter. Thus, for the partic- 
ular case under consideration, as shown in Fig. 5-6, the 
time of first ignition is dependent on items such as igniter 

mass flow rate (m. ) , ignition temperature (T . ) , 

xg ps f xg 

and convective heating rate, but the flame spreading rate is 

a motor characteristic which is largely independent of erosive 

burning, T . , and convective heating rate. Since the 
3 ps,ig 

insensitivities of flame spreading rate to erosive burning 

and T . (as shown on Fig. 5-6) are unexpected and not 
ps t xg 

fully explained, they should not be considered as generaliza- 
tions. 
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TIMESj 
0.014000-1 
0.06000-2 
0. 08000-3 
0. 10000-4 
0.50000-5 


V IGNITION TIME 



FIG. 5-3 STATIC TEMPERATURE ALONG PORT DURING IGNITION TRANSIENT 
(SRB1) 
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U. UU .1 a U U. - 't uu -x « u U 

01 S. FR0 M FORE-END OF PROP. CfUXlO 3 3 

FIG, 5-4 MACH NUMBERS ALONG PORT DURING IGNITION TRANSIENT 
(SRR1) 
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V IGNITION TIME 


TIMESs 
0.04000-t 
Q. 06000-2 
v 0. 08000-3 
U,L0000-4 
0.50000-^5- 



.00 0.80 1.60 2.140 3.20 14.00 

UIS. FROM FORE-END OF PROP. CfUXiO 3 } 

FIG. 5-5 STATIC PRESSURES ALONG PORT DURING IGNITION TRANSIENT 
(SRB1) 




FIG, 5-6 IGNITION TIMES ALONG PROPELLANT FOR A RANGE OF MOT 
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By studying the flow field development characteristics 
of Figs. 5-1 through 5-5, the reader can get a good under- 
standing of how flame spreading is implicitly an output of 
the p(x,t) model, i.e., 

- The hot gases from the igniter heat the propellant 
as they flow down the port. 

The rate at which the propellant is heated rapidly 
decreases in the direction of flow because igniter 
gases rapidly give up their heat and require a 
significant time to reach the nozzle end of the 
motor . 

- After the head end of the grain ignites, the flow 
rate of hot combustion gases along the port begins 
to accelerate and, thereby, accelerates the heating 
of the preheated (but unignited) propellant. 

As the flow rate increases, the acceleration of the 
combustion gases becomes one of the limiting factors 
and flame spreading rate becomes largely a character 
istic of the motor and not the igniter. 

As the hot combustion gases are driven down the 
port, the propellant is progressively heated to its 
ignition point, which is to say that flame spreading 
is described by successive ignitions. 

The intrinsic capability of the p(x,t) model to calculate 
flame spreading rates is one of its main features. 

Figure 5-7 is included to show a case in which the 
ignition delay is prominent. However, note that after 
first ignition occurs the p vs t characteristics are 
similar to those shown in Fig. 5-1, merely delayed by 0.1 
seconds. 


iEPRODUCBILipf OF 

PAGE IS POOR 



PRESSURE, ATMOSPHERES JX10 2 ) 

.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 
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FIG. 5-7 EFFF.CT OF 27% DECREASE IN IGNITER FLOW RATE AND 25% DECREASE 
IN CALCULATED CONVECTIVE HEAT TRANSFER COEFFICIENT, COMPARED 
TO CONDITIONS OF FIG. 5-1. 

(SRBA) 
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6.0 SOLID ROCKET BOOSTER CONSIDERED AS A SEGMENTED MOTOR 

The motor configuration considered in this section is 

intended to approximate that of the Solid Rocket Motor for 

25 

the Space Shuttle Booster, which (as shown in Fig. 6-la) 
consists of four segments. Since the burning surface area 
in the circumferential slots is inhibited, the burning 
surface associated with each circumferential slot is 
relatively small. The output associated with each slot is 
described in Table 6-1. 

As shown in Table 6-2 and Table 6-3 (parts 1 and 2) , 
the input is identical to that for the monolithic motors 
except that the geometry of the slots is specified as called 
for in Section 3'. Similarly, the output format as shown in 
Table 6-3 (parts 3 through 5; is the same as for a monolithic 
motor except that the state and contribution of the slots is 
included as described in Table 6-1. 

Figure 6-2 shows the pressure versus time at four 
positions along the grain. A direct comparison of the 
results for the segmented configuration (SRB/SEG1) with the 
results for the monolithic configuration reveals that the 
overall pressure response is very similar. However, as 
shown in Fig. 6-3 the flame spreading characteristics are 
influenced by the segments. These effects are attributed 
to the leading edge of each segment disturbing the boundary 
layer and thereby increasing the convective heat transfer 
coefficient. 

The small influence of the slots on the overall ignition 
transient is not surprising since the surface areas and 
volumes associated with the slots are relatively small 
fractions of the totals, for the particular design considered 

Figures 6-4, 6-5, and 6-6 show the pressure, gas 
velocity, and gas temperature of the three slots. Figure 
6-5 illustrates that the flows entering and exiting the 
slots are closely coupled to the pressure waves that occur 
in the main chamber. The flow oscillations in the slots are 




(a) Space Shuttle Solid Rocket Motor 


FWD SEGMENT 



(b) SRM Propellant Grain Design 


FIG. 6-1 TYPE OF SEGMENTED GEOMETRY CONSIDERED IN RUN SRB/SEG1 
(FIGURES TAKEN FROM AIAA PAPER NO. 75-1170, REF. 25). 
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Table 6-1 

OUTPUT ASSOCIATED WITH SLOTS 


Each time the output is printed, the following information 
is given for each slot: 


N 

n 

increment containing slot - 

PSLOT(N) 

p slot,n 

pressure in slot 

atm 

TSLOT (N) 

T 

slot,n 

temperature of slot 

K 

VSLOT(N) 

V . 

slot,n 

volume of slot 

3 

cm 

SWIDTH (N) 

w . , 
slot,n 

width of opening 

cm 

VELOCITY (N) 

v slot, n 

velocity of gases at 
slot/main stream 
interface. When > 0, 
flow is exiting. 

cm/sec 

MSLOT/MCH 

(pu) , ./(pu) . ratio of mass flux from 

slot to mass flux in 
main stream 


Example of 

output associated with each slot: 


PSLOT (ATM) 

TSLOT (K) 

VSLOT (CM3) SWIDTH (CM) VELOCITY (CM/SEC) 

1.05 

300.23 

0.49817E 06 5.8420 -1924 

.03 


Miscellaneous Output: 

A) Precedes first time dependent output: 

DELTAT=0 . 1667E-03 CSOUND=0 . 1061E 06 CPG=0.5872E 00 R=Q.3005E 04 

Initial time step Speed of sound at Specific heat 

seconds flame temperature of combustion 

conditions, cm/sec gases, cal/g-K 

B) Part of first time dependent output: 

PROPELLANT AREA, CMSQ=28 27989 . 00 CHAMBER VOLUME, CM** 3=75151920 . 0 

A b A dx 

Jo P 



reproducibility of the 
0RR5SNAL PAfil IS POOR 


51B/SEG1 2* 

tnu 

thax*.002 

TPRIKT=.001 
DELTAT*. 00016 
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TABLE 6-2 

INPUT DATA FOR SEGMENTED MOTOR (SRB/SEG1) 

BRO TPSCRI = 850 SAHPLE FOB IUSU 


HDELX*21 
tJHIT=-2. 
RP*PXT*0 
NPNPXT*2 
AT=2327. 

XP*3. 
IG*1341.3 
XB* 1377. . 

GAS A*1. 1363 
W =2 8. 21 
TFREF=3361. 
TIGN=2650 . 
RtJFS0R*0.01 
FRPR=0.00 1 1 
ROPR = 1.758 
CPR=. 3 
TPSCRI*750 . 
TPSCRT=85 0. 
RRBF-1.0783 
BRSXP - 0.35 
EBC*0. 

CH*.98 
DE=1 45.64 
ALFAD=12. 31 
RIGYAB-6 
RAPDVX • 22 


INPUT DIMENSIONS 
ARE IN INCHES 


NOTE: LEAVE COLUMN 

1 BLANK IN ALL 
NAMELIST INPUT 


BHD 

0. 

.020 

.250 

.340 

.450 

.50 

3. 


25. 

560. 

560. 

300. 

100 . 

0. 

2554. 



INPUTS IN ADDITION 
TO THOSE REQUIRED 
FOR MONOLITHIC MOTOR 


1018.92 1018.92 r 


133.85 

2554. 

1018.92 

1018.92 I 

199.3 

2870. 

608. 

608. 

244.7 

3187. 

202.83 

202.83 1 

330.1 

3167. 

202.63 

202.83 

395.6 

3187. 

202.83 

202.83 

461.0 

3187. 

202.83 

202.83 

526.4 

3187. 

202.83 

202.83 

591.9 

3187. 

202.8 3 

202.83 

657.3 

3187. 

202.83 

202.83 

-722.7 

3187. 

202.8 3 

202.83 

788.1 

3187. 

202.93 

202.83 

853.6 

3187. 

202.83 

202.83 

919. 

3187. 

202.83 

202.83 

584.4 

3153.7 

202.83 

202.83 

1849.9 

3127. 

202.83 

202.83 

9115.3 

3203. 

202.8 3 

202.83 

1180.7 

3278. 

202.83 

202.83 

1246.1 

3354. 

202.83 

202.83 

1311.6 

3709. 

202.83 

202.83 



202.83 •' 

'■ 202.63- 


Blot 

XSLOT 

361. -'f 


681. 


1001. 


®S,Blot 

DELABS 


<$V 


slot 

DELVS 


411, 34448. 


2660 . 266H. 


3212. 


"■lot 
WSLOT | 

S*» 


2.3 


A fcl 
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TABLE 6-2 (CONTINUED) 

INPUT DATA FOR SEGMENTED MOTOR (SRB/SEG1) 


SHAME 

!7, Max=. 01 
TPBINT=0. 002 
DEL ?AC=1 . 3 

SEND 

SNAMF 

TPRtNT=.005 
TMAX-. 1 1 
DEL FAC=2 * 

SEND 

BNAHE 

?M A X= . 2 

SEND 
SNA ME 

TM A X= . 4 

TPOTJpr, 0 1 

DEL r AC=1 . 5 

6 EN D 
SNAME 

Ti*AX=.6 
A X- . 5 
DEL FAC=2 . 

SEND 



HIGH PERFORMANCE ROCKET MOTOR IGNITION TRANSIENT PREDICTION PROGRAM - SEGMENTED - FEB 1976 
INCLUDES: SPATIAL 6 TIME DEVELOPMENT OF P.U.ET AND FLAME SPREADING 

SRB/SEGl 21 NO ERO 1PSCR I = t}60 NEW FCT2 

GENERALI LIMITS. CODE WORDS, 6 INITIAL CONDITIONS: 

TMAX LAMBDA TPRINT P2DNE DECT AT NDLLX 

'- #.002 1. 00 0.001 7777.00 0.0001600 21 

TP I PAM UNIT NINER1 NPNPXT 

298.0 1.00 -2. 0 2 

MOTOR CONFIGURATION PARAMETERS: 

ATT XP XG XL 

2327.00 3.00 1341.30 1377.00 

COMBUSTION GAS PROPERTIES 
GAMA W TFREF TIGN 

1.136 28.21 3361.0 2660.0 

FRICTION TRANSFER AND HEAT TRANSFER CQCF.: 

RUFSUR DDRG ODHC 

0.0100 1.00 1.00 


i»fffiRELLANT PROPERTIES: 




FJWI 

ROPR 

CPR 

TOREF 

SIGP 

TPSCR1 

0.001 iO 

1.758 

0.30 

296.00 

0.00200 

660.00 

BURNING 

RATE RELATIONSHIP 

“ 



RREF 

PREF 

BREXP 

EBC 

EbEX 


-1.078 . 

08.08 

0.3500 

0.0 

126.000 


PARAMETERS USED 

IN THRUST 

CALCULAT IONS: 


DE 

CM 

ALFAD 

EROAT 

EROEXP 


1*5.6* 

0.980 

12.31 

0.0 

0.60 



JOMTTER: MASS DISCHARGE VS. TIME TABLE: 
TIME '■■■'■ FLOW RATE 
0.0 26.0 

0.0200 660.0 

0.2500 560.0 

- ' 30.3*00 300.0 

: ~ 0.460 0 100.0 

SE.5SW p 0.0 



TABLE 6-3 (PART 1) 

OUTPUT FROM SEGMENTED MOTOR RUN (SRB/SEG1) 
- INPUT - 



PROPELLANT GEOMETRY TABLE: 


LENGTH PORT AREA 

PROPELLANT 

WETTED 

UPSTRM OF 

BURN AREA 

VOLUME 

WIDTH 



PERIMETER 

PERIMETER 

SLOT 

OF SLOT 

OP SLOT 


l Tj'ii’ l ^ r • 

2554 .OO 

1018.92 

1018.92 

0.0 

0.0 

O. . 

0.0 ••• =\ . 

4**43 

2554.00 

1018.92 

1618.92 

0.0 

0.0 

e. '■ '• 

. o. e ===--== 

133.86 

2554.00 

1018.92 

1018.92 

0.0 

0.0 

0. 

0.0 

199.30 

2870.00 

608.00 

608.00 

0.0 

0.0 

0. 

0.0 

264.70 

3187.00 

202.83 

202.83 

0.0 

0.0 

o. 

0.0 

3311.10 

3487.00 

202.63 

202.83 

361.00 

411.0 

30400. 1 

2.30 

»&.u 

3187.00 

202.83 

202.83 

0.0 

0.0 

o. 

0.0 / • ' 

4612W 

3187.00 

202.83 

202.83 

0.0 

0.0 

0. 

0.0 . . ' > 

526. AO 

3187.00 

202.83 

202.83 

0.0 

0.0 

o. 

o.o 

591.90 

3187.00 

202.83 

202.83 

0.0 

0.0 

0. 

0.0 

657.30 

3187.00 

202.83 

202.83 

681.00 

2660.0 

28600. 

2.30 

: '' ' T«*TO - v " 

3107.00 

202.83 

202.63 

0.0 

0.0 

o. 

o.o 

v TSS.IO 

3187.00 

202.83 

202.63 

0.0 

0.0 

0. 

0.0 . • • . 

as 3. fco 

3187.00 

202.83 

202.83 

0.0 

0.0 

0. 

0 . 0 .,. : .. 

919.00 

3187.00 

202.83 

202.83 

0.0 

0.0 

0. 

0.0 

984.40 

3153.70 

202.83 

202.83 

1001.00 

3212.0 

28600. 

2.30 

1049.90 

3127.00 

202.83 

202.63 

0.0 

0.0 

0. 

0.0 

■;i. »T*230 

3203.00 

202. 83 

202.83 

0.0 

0.0 

O. 

0.0 “ • 

% 1180.70 

3278.00 

202.83 

202.83 

0.0 

0.0 

O. 

0*0 

1246. LO 

3354.00 

202.83 

202.03 

0.0 

0.0 

0. 

0*0 

131 1 .60 

3709.00 

202.83 

202.83 

0.0 

0.0 

0. 

0.0 

1377.00 

4500.00 

202.83 

202.83 

o.o 

0.0 

0. 

0.0 


BURNING RATE TABLE - R VS P 


MDATA NPPR 

0 2 

PRESSURE RATE 

***4GEOtt£TRY CONVERTED FROM INCHES TO CENTIMETERS 


length 

PORT AREA 

PROPELLANT 

PERIMETER 

WETTED 
PER I ME TER 

UPSTRM OF 
SECT 

BURN AREA 
OF SLOT 

VOLUME 
OF SLOT 

WIDTH 

OF SLOT 

7.62 

16477.38 

2588.06 

2588 . 06 

0.0 

0.0 

0. 

0.0 

173.81 

16477.38 

2588.06 

2588.06 

0.0 

0.0 

0. 

0.0 

340.00 

16477.38 

2588.06 

2588.06 

0.0 

0.0 

O. 

0.0 

506.22 

18516.08 

1544.32 

1544.32 

0.0 

0.0 

0. 

0.0 

672.34 

20561 .24 

515.19 

515.19 

0.0 

0.0 

0. 

0.0 

838.45 

20561.24 

515. 19 

515.19 

916.94 

2651.6 498166. 

5.84 

1004.82 

20561 .24 

515.19 

515.19 

0.0 

0.0 

O. 

0.0 


TABLE 6-3 (PART 2) 

OUTPUT FROM SEGMENTED MOTOR RUN (SRB/SEGl) 
- INPUT GEOMETRY - 


TABLE 6-3 (PART 3) 

OUTPUT FROM SEGMENTED MOTOR RUN (SRB/SEG1) 

- DURING FLAME SPREADING - 

TIME =0*9000E-01 I 208 DEETATl DEL. TAT *0.5000E-03 


h DISTANCE PRESSURE T MACH RATE BR/BRN TAUB FR1CT. CUNY. M/A tSURF PORT 



CM 

IN 

ATM 

PS 1 A 

DEG K. 

uu. 

cm/sec 


CM FACTOR COEF. G/SEC CMSO 

K 

CMSO 



? 
















7 -to 

3-0 

15.83 

232»b 

2t»t>3*3 

0.078 

0.65 

1 . 000 

0.0208 

0.0 

0.0315 

15.16 

850. 

16531.1 


1 

iT3. a 

68.4 

15.09 

22 t . 7 

2-938.8 

0.233 

C *64 

1 .00 0 

0.0160 

0.0 

0.0645 

40.87 

850. 

16518.8 


2 

340.0 

133.9 

13.40 

i9t>.9 

3C74.2 

0.415 

0.61 

1 .000 

0.0130 

0.0 

0.0880 

63.06 

850. 

16511.1 


3 

SOto.2 

199.3 

1 0.64 

1 So . 3 

3u64. 1 

it . 623 

0.56 

1 .000 

0.0100 

0.0 

0.0920 

75.07 

850. 

18530.9 



264.7 

9.01 

132.3 

3066.8 

0.749 

0*53 

1 .ooo 

0.0076 

0.0 

0.0808 

7to.Sto 

850. 

20565. 1 


*> 

838 -to 

33c- 1 

8.81 

129.4 

31 05.6 

0.759 

0.53 

1.000 

0.0068 

G.O 

0.0785 

75.45 

850. 

20564.7 



1004 . a 

395.6 

t4 ■ btj 

125 <• b 

3143.2 

0.774 
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TABLE 6-3 (PART 4) 
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(PART 5) 
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FIG, 6-2 PRESSURE DEVELOPMENT ALONG PORT OF SEGMENTED MOTOR 
(SRR/SEG1 FOUR SEGMENT MOTOR) 
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FIG, 6-3 FLAME SPREADING TIMES OF SEGMENTED MOTOR COMPARED WITH 
FLAME SPREADING TIMES OF MONOLITHIC MOTOR 


(SRB/SEG1 COMPARED TO SRB1) 
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Fig. 6-5 Following ignition of burninq surface in slots, 

the gas velocity entering and exiting the slot is 
very sensitive to pressure oscillation in the 
main chamber. 
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Fig. 6-6 Mean gas temperatures in slots lag the 
gas temperatures in the main stream. 
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complex, since each slot has its own characteristic filling 
and venting time. During the flame spreading interval, the 
net flow is entering the slots. After ignition of the burn- 
ing surface area in the slots and during the early stages of 
chamber pressurization, the small amount of mass generated 
by the relatively small burning surface in each slot serves 
to partially pressurize the volume in the slot. The two 
downstream slots with their larger surface areas produce a 
net outflow of gas during the early phase of chamber pres- 
surization, whereas the upstream slot begins to produce an 
outflow toward the end of chamber pressurization. As shown 
by the temperature differences (between the chamber and 
slot gases) in Fig. 6-6, the time required to flush the 
warm gases from the slots is appreciable. 

A comparison of the pressure oscillations predicted for 
the monolithic motor (Fig. 5-1) with those of the segmented 
motor reveals that the volumes between the segments and mass 
flow entering and exiting the slots serve to damp the longi- 
tudinal waves. This illustrates a very important intrinsic 
feature of the solution — the ability to analyze certain 
aspects of longitudinal combustion instability in segmented 
motors. Indeed, a logical follow-up of the present study 
is to extend the computer program and use it as part of 
a study to determine the extent to which the slots in the 
present SRB design reduce the possibility of axial mode 
combustion instability. 

As shown in Fig. 6-1, the nozzle of the head-end igniter 

is very close to star-point tips in the first segment. It 

2 6 

has been pointed out that treating the head-end region as 
a uniformly heated port (as in the present analysis) results 
in unrealistically long induction periods (see Fig. 2-1 
definition) . A better simulation would be obtained by con- 
sidering at least three zones of heating and ignition in the 
first segment: (1) the more intense heating and rapid igni- 

tion of the star-point tips, (2) the increasing heating 
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rate in the axial slots as the gas generated by the burning 
star-point tips augments the axial flow of hot gases, and 
(3) the conventional heating of the aft portion of the first 
segment. Thus, the net effect will be to reduce the pre- 
dicted induction time which is a major contributor to the 
uncertainty associated with motor-to-motor variation in 
time-to-full-chamber-pressure. Recall that the results of 
Section 5 showed that flame spreading rate is relatively 
insensitive to the variations that may occur within the motor, 
whereas time to first ignition is affected greatly by small 
variations . 

The analysis of the interaction of the slot flow and 
the main chamber flow is based on the assumption that the 
slots have been properly designed and that the aerodynamic 
interactions between the two flows will not produce appre- 
ciable pressure drops across the slots. Reference 27 
describes analytical and experimental methods for deter- 
mining the severity of aerodynamic interactions and pre- 
scribes a design approach to eliminate the undesirable 
affects. The user of the program is cautioned to insure 
that the particular slots being analyzed conform to the 
assumptions and conditions stated in Section 2. 
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7.0 CONCLUSION 

The analytical model described in this report provides 
for the first time a means of analyzing the complexities 
of ignition transients and pressure peaks of large, high- 
performance, segmented SRB's. This is accomplished by 
accounting for (1) the temporal and spatial development of 
the flow field set up by the igniter discharge, (2) 
ignition and flame spreading coupled to chamber flow, (3) 
the large velocity, pressure, and temperature gradients that 
occur during the early phases of ignition, and (4) the 
interactions that combine to produce peak pressures (i.e., 
compression of chamber gases during pressurization, erosive 
burning, and mass-added effect of igniter discharge) . 

As the pursuit of higher performance rocket motors 
continues and the present configurations are upgraded 
(i.e., by increasing loading density, using high-performance 
propellant, and extending ambient temperature range) , the 
methods described in the report will enable analysts to 
predict maximum pressures with much greater confidence 
than was heretofore possible. 

The calculated results in this report are intended to 
introduce the user to the computer program and are not 
intended to be a comprehensive investigation of a particular 
SRB design. Nevertheless, the results of Section 5 and 6 
produced several interesting insights. Once the head-end 
section of the SRB ignites, the subsequent flame spreading 
rate is a motor property that is largely unaffected by 
changes in the igniter, the erosive burning contribution, 
and the propellant ignition temperature; whereas, the time 
of first ignition is strongly dependent on igniter charac- 
teristics and propellant ignition temperature. The contri- 
butions of the burning surface areas and volumes in the 
slots of the present SRB design do not greatly influence 
the times of the important ignition events. The interactions 
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of the main chamber flow with the slots tend to damp axial 
pressure waves which are set up toward the end of flame 
spreading. 

The ignition prediction and analysis procedures 
described in the report can be applied to a wide variety 
of rocket motor applications and are not limited to large, 
segmented boosters. Special purpose applications can be 
achieved by extensions to the present program, e.g., aft- 
end ignition; flame spreading delays produced by restrictors 
(e.g., used to reduce ignition spikes); ignition of degraded 
or wet propellant surfaces; case expansion (and A /A 

p t 

increases) during pressurization of low safety factor 
motor cases; and nozzle closures. In its present form the 
program contains many of the elements necessary for the 
analysis of axial mode combustion instability in segmented 


motors . 
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APPENDIX A 


Performance Parameters 


The motor performance equations are presented to specify 
the manner in which the motor parameters are applied. A 
compatible set of combustion gas properties is established by 
starting with calculated combustion gas properties (obtained 
from thermochemistry computer programs) and then adjusting 
the calculated properties to obtain agreement with measured 
mass flow rates and thrusts. Careful attention should be 
given to accounting for items such as heat loss, combustion 
efficiencies, and ambient temperature. The measured value 
of c* that is normally used to calculate the mass discharge 
is established from motor p vs t data and weight of pro- 
pellant burned, 


w pr , by the following equation 


9K 


E, stag 


raeas 


w 


(A-l) 


pr 


The quantity c* is used to relate instantaneous mass 
in© ci s 

flow rate, m , to throat area, A fc , and chamber pressure. 


^E, stag 


m n ^t^E , stag^^meas 


(A- 2) 


From isentropic one-dimensional compressible flow theory, 
mass discharge rate can be expressed in terms of the gas 
properties as follows. 


Y 


1/2 


Y + 1 


Y + 1 

2(y - IT 


m P E,stag A t g 




gR g T E,stag /W 


(A- 3) 
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Ideal ly, c* is defined in terms of the chamber gas proper- 
ties as 



A 


gR t £ /w 


X/2 


2 

y~TT 


Y + 1 

2Ty -1) 


J 


(A- 4 ) 


i.e., the chamber temperature is the same as the flame 
temperature under no loss, equilibrium conditions. 

To have a compatible set of relationships, the consist- 
ency between c* __ and T , established by Eq. (A-5) 

IuSuS Cil / aC u 

is maintained throughout the analysis.* 


T 

ch,act 


W 

gR, 


C* Y 

meas ' 


1/2 


I ± 1 

277 - 1 ) 


.Y + lj 


(A-5) 


Thus, the input value of T^ used in the program should 

be T ch,act* 

The ratio of specific heats, y, is an average effect- 
ive value for the motor and nozzle. When possible, this value 
was determined from experimental data. Because of shifting 
equilibrium during flow through the nozzle, this value is 
little more than empirical correlation constant. When sat- 
isfactory experimental data is not available, an effective 
value based on thermochemical calculations can be deduced 


♦Using T ^ act has the effect of lumping other losses (i.e., 
nozzle discharge* and two phase flow) with combustion efficiency 
losses. A more thorough approach to establishing T ^ t could 
be based on the Standard Performance Prediction program recently 
developed by the U.S.A.F. Rocket Propulsion Laboratory or the 
Thiokol Performance Prediction program recently developed by 
the Wasatch Division of the Thiokol Corporation. 
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from the relationship between theoretical gas properties 


at chamber and exit condictions. 



Y - .1 

T . 
ch „ 

[ p ch] 

Y 

T 

ex 

p exj 



or between the theoretical gas properties at chamber and 
nozzle throat conditions. 



f 



(A— 7 ) 


Thrust is expressed in terms of a thrust coefficient as 


F C FAm P E, stag A t 


(A- 8) 


The particular thrust coefficient used in the program is 


C FXm C m 


C FV X N + (1 “ Vp 
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ex 


E,stag 


ex 


P. 
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P ^ex <*- 9 > 

P E, stag 


where 


C FV + 


Y + 1 


2y: 


Y " 1 


lY + 1 J 


Y " 1 
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Y - ll 
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P ex 
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. E, stagj 

• 



1/2 


+ e 


ex 


eX P E, stag 


(A-10) 


and 


X N = (1 + cos « n )/2 


(A-ll) 


The value of motor coefficient, C m , used in the equation 
is for nozzle-end stagnation pressure and should not be con- 
fused with values that were determined from the relationship 
between head-end pressure and measured thrust. Also, note 
that does not multiply the last term in Eq. (A-9) . 
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APPENDIX B 

MODIFIED* C.G.S UNITS FOR PROPELLANT COMBUSTION CALCULATIONS 

and 

APPROXIMATIONS TO GAS PROPERTIES 


Gas Law: 


P/P = (R /W)T 
p - [atm] 
p - [g/cm 3 ] 


T -[K] 

W -[g/g-mole] 


(B-l ) 


R = 82.05 

g 


cm atm 


Lg-mol-KJ 


Speed of sound: 


c = (gyp/p)^^^ = (gyTR^/W) 3 ’ // ' 2 [cm/sec] 


1.013 x 10' 


(B-2) 


_atm-cm-sec 

Viscosity: Empirical fit for typical combustion products 

.-6 rT 0 . 5„,0 . 6 


p = 1.185 x 10 W ’ T * [g/cm-sec] 

Thermal conductivity of low density combustion gases: 
Polyatomic Eucken equation: 


(B— 3 ) 


X = 


r R 

c + - 
[ C p 4 W J 


\i [cal/cm-sec-K] 


(B-4 ) 


c p - [cal/ g-K] 


R = 1.987 
u 


cal 

JOl-K 


g-mol 


Mass flow rates: 


* _ P E , stag^t^ 

rp c A b 

r -[cm/sec] , 


A fa - [cm ] 


[g/sec] 

p -[g/cm 3 ] 
p - [atm] 


(B-5) 


♦Pressure in atmosphere rather than dynes/cm . 
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(gTR /W) 1//2 

Y — units same as speed of sound (B-6) 


Reynolds number: 




(B— 7 ) 


p = 


T (R^/W) 
y 


Specific heat: 


R 

c « _JL_ -H 
p,g y-i w 


(B-8) 


(B-9 ) 


Internal energy: 


u 


e + 2j~ Ical//g] 


(B-10) 


Prandtl number (Eucken formula) : 


c u 

Pr = 1 = J?1 = 
a 


X “ c +1.25R u /W 


(B-ll ) 
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